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ADDITIONAL NOMENCLATURE 


See Reference 1 for the basic nomenclature. An asterisk before the 
symbol indicates an input variable. All input subscripted and non- 
s ibscripted variables are listed separately. 


English 

Symbol Definition Units Used 


A GS1’ A GS2’ A GS 

Cross-section area of a wagon-wheel grain 
associated with the first, second and both 
sets of grain points, respectively. 

in 2 

#c* 

ig 

Characteristic exhaust velocity of the igniter 
propellant 

ft/sec 

#h *h 
1’ 2 

The half width of the first and second set of 
points, respectively, of a wagon-wheel grain 

in. 

*1 

go 

Integer designating option of ignition transient 
calculation. 


*1 

po 

Integer designating option of plotting results 


*1 

wo 

Integer designating option of inert weight 
calculations 


* K a ,#K b 

Empirical constants in the characteristic 
velocity versus chamber pressure linear 
relations (Eq. B^c) used for ignition transient 
calculations 

ft/sec , 
ft/sec-psia 


The length of the two parallel sides of the 
first and second set of points, respectively, 
of a wagon-wheel grain 

in. 

*N(jj) 

Integer designating whether or not a specific 
output plot is desired 


#p 

big 

Value of chamber pressure P c at which a 
noz2le blowout plug in the main motor is 
ejected 

psia 
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English 

Symbol 

Definition 

Units Used 

*P . 
mig 

Maximum chamber pressure attained in the 
igniter 

psia 

h 3 

Distance from center of curvature of a 
spherical end of circular perforated grain 
to the burning surface associated with 6^ 

in. 

*R. 

IS 

Average regression rate of the first half 
of the igniter pressure-time trace 

psia/sec 

S G1 ,S G2 ,S G 

Perimeter of a wagon-wheel grain associated 
with the first, second and both sets of 
grain points, respectively. 

in. 

'**11 ’**12 

Time for the igniter, pressure to reach 
maximum value and to’ decay to 10 percent 
of its maximum value, respectively 

sec. 

* U fs 

Flame spreading speed during ignition 

in/sec 

* V ciT 

Initial volume of chamber gases associated 
with tabular input 

in 3 

X 

o 

V* 

*< 

o 

Coordinates of tips of wagon-wheel grain 

in. 

V’V 

Greek 

Symbol 

Coordinates of intersections of burning 
slant sides of wagon-wheel grain points with 
fillet arcs (origins at fillet centers) 

in. 

#(V 

V 2 

The angle between the slant sides of a wagon- 
wheel grain point and the center line of the 
point for the first and second set of points , 
respectively 

radians 

*At. 

ig 

Time increment for ignition transients 

sec 

6 fw 

Angular location of fillet centers with 
respect to radial centerline of wagon- 
wheel grain points 

radians 

X 

Volumetric loading density; i.e., initial volume 
occupied by propellant divided by empty case 
volume 
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Greek 

Symbol 

x 

ww 

Subscript 

ig 


Definition 


Units Used 


Web thickness of wagon-wheel grain. 


in. 


ignition 
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I. INTRODUCTION AND SUMMARY 


This report presents the results of research performed at Auburn 
University during the period from December 15, 1972, to April 15, 1972, 
under Modifications Nos. 8, 9 and 10 to the Cooperative Agreement, dated 
October 6, 1969, between the NASA Marshall Space Flight Center and Auburn 
University. The research is an extension of that accomplished under 
Modification No. 6 to the same agreement which culminated in the develop- 
ment of a simplified computer program for preliminary design and performance 
analysis of solid-propellant rocket motors (SRMs) as reported in Reference 1. 

The extension adds the following capabilities to the simplified program 
as program options . 

1. Treatment of "wagon-wheel" cross-sectional propellant 
configurations alone or in combinations with circular 
perforated configurations . 

2. Calculation of ignition transients with the igniter treated 
as a small rocket motor. 

3. Accurate representation of spherical circular perforated 
grain ends as an alternative to the conical end surface 
approximation used in the original program. 

4. Graphical presentation of program results using the 
"CalComp 663" digital plotter. 

In addition, inert weight calculations are made optional in the new 
program whereas they were required calculations in the original program 
necessitating specification of input values. Also, the initial and final 
gaseous chamber volumes are eliminated as input requirements in the new 
program. These values are now calculated from the other program inputs. 
Finally, the computer program has been modified in detail for the purpose 
of simplifying the program and minimizing the computation time. The net 
result of the extensions, however, is an increase in compilation time to 
a maximum of about 2 minutes and 20 seconds and an increase in execution 
time to about 20 seconds in the IBM 360 computer. Almost all of the 
increase is attributable to the ignition transient calculations . The new 
program uses approximately 19 ,000 words on the IBM 360 . 

The present report is prepared as an addendum to the original report. 
Familiarity of the reader with Reference 1 is assumed. Thus, the extensive 
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notation of the original report is not repeated; only the nev notation is 
given. The nev input variables or deletions are discussed in Section II. 

The organization of the entire report follows that of the original report 
and equations are numbered consecutively with respect to the original report 
to provide for insertion at the appropriate locations in the original 
analysis. Figures and tables are numbered similarly according to whether 
they replace or supplement those of the basic report. 

Flowcharts of the two new subroutines (ignition and graphical representa- 
tion) are presented in Section IV but no attempt has been made to modify the 
flowcharts of the main program because the logic of the modifications is a 
straightforward extension of that of the original program and is made clear 
from the mathematical analyses presented in Section III. The new data card 
format is shown in Table IV-1 and the complete extended program is listed in 
Table IV-2. Program statements that must be removed in order to delete the 
CalComp plotter compilation requirements are indicated by check marks (►) in 
Table IV-2. Removal of these statements is necessary if the computer is not 
equipped for "CalComp" plotting. However, if alternate plotters are avail- 
able, generally only the plotting subroutine need be replaced. 

Test cases are given in Section V to illustrate (l) ignition transient 
calculations, (2) spherical closure effects with circular perforated grains, 
(3) graphical representation of program result, and (4) computation for 
"wagon-wheel" cross-sectional grain configuration. 



II. DISCUSSION OF INPUT AND OUTPUT 


In this section each new input variable is defined in the order in 
which it is encountered in the program. Deletions or other changes in 
original variables are similarly treated. As in Reference 1, appropriate 
additional discussion of the variable and typical numerical values are 
given. New output variables of the program are also identified and defined. 


Users Options 

I (IGO) 0 

S° 

1 

I (IWO) 0 

wo 

1 

I (IPO) 0 

po 

1 

2 

3 

N( j j ) ( numplt(jj )) 

0 

1 

1 

2 

3 


For no ignition transient computations 
For ignition transient computations 
For no inert weight computations 
For inert weight computations 
For no plots 

For plots of equilibrium burning only 

For plots of ignition transients only 

For plots of both ignition transients and equilibrium 
burning. 

An integer designating whether or not a specific 
output plot is desired: 

If a specific plot is not desired 

If a specific plot is desired 

The order of specification of NUMPLT(JJ) is 
as follows : 

PHEAD versus T 

PONOZ versus T 

PHEAD and PONOZ versus T (superimposed) 
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RHEAD versus T 

5 RNOZ versus T 

6 RHEAD and RNOZ versus T (superimposed) 

T SUMAB versus T 

8 SG versus T 

9 SUMAB and SG versus T (superimposed) 

10 F versus T 

11 FVAC versus T 

12 F and FVAC versus T (superimposed) 

13 VC versus T 

14 SUMAB versus Y 

15 SG versus Y 

1 6 SUMAB and SG versus Y (superimposed) 

Primary Basic Motor Dimensions 

V . ,V (VCI,VCF) Initial and- final volume of chamber gases, respectively. 
'~ 1 These are deleted as input variables in the new program 

except for volume associated with tabular input. 

See V ciT below. 

Ignition Characteristics . Not required if ignition calculations are not 

desired (IGO = 0) 

K ,KL (KAjKB) Empirical constants in the characteristic velocity versus 

a chamber pressure relationship for the main motor (ft/sec, 

ft/sec-psia) . A linear relationship is as sinned. Data 
are obtained from thermochemical calculations. 

Flame spreading speed (in/sec). Reference 2 cites a 
value of 3650 in/sec for one composite type propellant 
with ammonium perchlorate oxidizer and aluminum 
additives. Values will differ for various specific 
propellants . 


U fs (UFS) 
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C* (CSIG) 
lg 


Characteristic velocity of the igniter propellant 
(ft/sec). 


P . 

mig 


(PMIG) 


t n (m) 


t TO (TI2) 


R. (RRIG) 

"*■ O 


At. (DELTIG) 

lg 



(PRIG) 


Maximum chamber pressure of igniter when discharging 
into the atmosphere (psia). This pressure is 
usually established based upon minimization of the 
igniter weight. It will ordinarily range from 80 
to 150 percent of the maximum operating pressure of 
the main motor. 

Time for the igniter pressure to reach its maximum 
pressure (sec). Values range from 0.010 sec for a 
small igniter to 0.050 sec for a large igniter. 

Time for the igniter pressure to decay to 10 percent 
of its maximum value when discharging into the atmosphere 
(psia/sec). Typical values range from 0.2 to 0.5 
seconds in smaller igniters to 0.5 to 1.2 in very large 
igniters . 

Average regression rate of the first half of the igniter 
pressure time trace (psia/sec). Most igniter traces 
are generally regressive or two-level with the higher 
level occurring first. The analysis considers only 
the first half of the pressure time trace, because 
the contribution of the mass flow rate of the igniter 
to total mass flow rate is insignificant beyond this 
region during a normal ignition. 

Time increment for calculation of ignition transients 
(sec). A value of 0.001 seconds is recommended for 
small motors and 0.005 seconds for very large motors. 

Value of chamber pressure P c at which a nozzle blowout 
plug in the main motor is ejected (psia). If no 
blowout plug is used, PBIG should be set equal to an 
estimate of the atmospheric pressure when ignition 
is initiated. 


Basic Properties of Weight Calculations 

These inputs as listed in the original program are no longer required 
if calculation of inert weight is not desired (IWO = 0) 
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Input to Establish the Program and Basic Grain Configuration 
and Arrangement 

Add: 3 for wagon-wheel (See Figure II-4) 

For extreme ends of a circular perforated grain only: 

0 If "both ends are conical or flat. 

1 If head end is conical or flat and aft end 
is spherical. 

2 If "both ends are spherical. 

3 If head end is spherical and aft end is conical 
or flat . 

If the end of a circular perforated grain is taken spherical, the correspond- 
ing angle THETACH or THETACH is not used in the program. 

Tabular Burning Surface and Port Areas (Not required for INPUT=2) 

V c i<r Initial volume of chamber gases associated with tabular 

input (in 3 ) 

Basic Geometry for Star Grains 

In order to simplify the conversion from the original to the extended program, 
the wagon-wheel is considered a type of star grain and each variable in this 
group as listed in the original program must be specified for wagon-wheel as 
well as standard star or truncated star grains. 

Special Geometry for Wagon-Wheel Grain (Not required for standard or 

truncated star grains (STAB = 1 or 2)) 

The web thickness of wagon-wheel grain (in.). The 
discussion under TAUWS for standard star grains is 
applicable. 

The lengths of the pairs of parallel sides of the first 
and second set of grain points, respectively (in.) 

(See Figure II-4). 

The angles between the slant sides and the center lines 
of the points of the first and second set of grain 
points, respectively (radians). The angles should 
not exceed tt/2 radians , 

The half-width of the star points (in.). HW must not 
exceed TAUWW (HW < TAUWW). 


T 

ww 

(LI ,1.2) 

a i ,a 2 

(ALPHA1 ,ALPKA2) 

h (Htf) 
w 


S (STAR) 
op 

C (COP) 
op 
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Additional Program Outputs 

The outputs listed below are in addition to those of the original 
program. However, those outputs related to the ignition transient calcula^ 
tions are only obtained when IGO =1. It should be noted that values of 
pressure (PONOZ and FHEAD) are based upon the assumption of choked flow 
through the nozzle and values of delivered thrust (F) and related values 
(e.g., CF) are based upon the additional assumption of the nozzle flowing 
full (no separation) . This introduces errors in the output parameters 
mentioned during the initial phase of the ignition. 


Outputs as Functions of Time 

t. (TIG) Operating time during ignition (sec,). This is calculated 

from the time of initiation of flame spreading. 

P c (i )(PCIG) Chamber pressure within the igniter (psia). This is 

s based on a piecewise linear pressure versus time relation 

before and after peak igniter pressure is encountered 
until the main motor pressure feeds back and controls 
the igniter pressure . 


Non-Time Varying Quantities 


A*. (ASIG) 
ig 

”W) (MIGAV) 

W ig.Ctot)' WIGT0T) 


Required igniter throat cross-sectional area (in 2 ). 

Average mass flow rate of the igniter discharge over 
the first half of the igniter burning time excluding 
time prior to Til (slugs/sec). 

Estimated total weight of required igniter propellant 
(lbs.). This is based on the integral of m-^ over the 
time interval from Til to T12/2 plus a somevfiat arbitrary 
allowance of two-thirds the stated integral to provide 
sustained heat transfer feedback during the final 
ignition phase and propellant for the pressure buildup 
phase within the igniter (see Eq. B9). 


V ci > V cf (VCI ’ 7CF) 


Initial and final volume of chamber gases , respectively 
(in 3 ) 


\ (LAMBDA) Volumetric loading density; i.e., initial volume 

occupied by the propellant divided by the empty case 
volume Cl). Note discussion after Eq. 37d. 



Ill . ANALYSIS 


In this section the extensions to the mathematical analysis of the SRM 
design problem are presented in a step by step procedure following the method 
of the original report (Reference l). The steps in the procedure are 
numbered consecutively with respect to those steps of the original report 
to permit insertion at the appropriate locations in the original analysis. 

In the few instances where one of the original analysis steps has been 
revised, the new analysis step bears the same number as the original step. 


Burning Surfaces and Port Areas 

End Geometry of Circular Perforated Grains 

AO.h. Select an option for the circular perforated grain end configuration. 

If both ends are conical or flat, set C = 0. 

op 

If head end is conical or flat and aft end spherical, set C = 1. 

If head and aft ends are both spherical, set C = 2. 

If head end is hemispherical and aft end conical or flat, C ^ =3. 


AU.b. (revised) For the perforations. 

If 2 + D. > D , set A = A, = 0, 
y i o’ Dnc Dpc 

If not and if <_ 5° and 

if C = 0 , compute 
op 

L_ = L„ . - y (cot & + 6 , ) 

Gc Gci cn ch 

if C =1 compute 
OP 

L = L_ . - { [D 2 - (D. + AD.) 2 ] 3 
Gc Gci o i l 


Go to step A 5 


- (D q 2 - (^ + AD ± + 2y) 2 ]^}/2 - y cot 0 


ch 


if C = 2, compute 


^ ■ ^Gci - - 'V - (D i * * )2 * 
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if = 3, compute 


V - Vi - UV-Vl 

- [D 2 - (D. + 2y) 2 ]**}/2 - y cot fl 


*bpc = " (D i + 2y,[1 0c - T t =S ! fOT Vo > 


Vo * ° f0r Vo -“'2 


If 6„ > 5° and 
G 

if C = 0 or 1, compute 
op 


Vo ’ ' (D i + ^‘Vi - y cot V 

- T £ - Is + tan(e G /2)]y} for A bpc > 0 


if C = 2 or 3, compute 
op 

Vo - * (D i + * ) <Vi - ■ttV-V’’ 5 

- ID 0 2 - (D i + 2y) 2 ]^}/2 - T a 

- IS + tan (0 G /2)]y> for > 0 

A, - 0 for A. <0. 

Dpc Dpc — 

A^.c. (revised) For the nozzle end surface. 

If e G < 5°, set = 0 

If 6. > 5° and C = 1 or 2 
G op 

R 3 = I(Di + ADi)/2 + L Gni sin 0^ cos © G 
- (sin 6_) ' {(D /2) 2 - I (Li + ADi)/2 

G O 

+ Vi sin 9 g ]2))s 
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Vc = * tL Oni - R 3 - y tan < e G /2)] + ADi >/ 2 

+ y + (sin e G )[(D o /2) 2 - (B 3 + y) 2 ]^ 

+ (R + y) cos 0 } for A >0 
i G one 

A. " 0 for A, <0 
one one 

If 6 n > 5° and C = 0 or 3 
G op 

Vn * ’ tL Gni ‘ y Cot (6 G + 6 <m ) - y tan (6 0 /2)][l>i+4Di 

* y * L Gni 3in 6 G + y 0SC (8 G * "ax' 1 6ln 6 on ] for Vc ’ 0 

A, =0 for A, < 0 
Dne one — 

Note: For 0^ > 5° » the entire "burning end. surface is treated as 

conical although a non-coni cal. surface area evolves when 

0 _ + 0 > ir/2 for C = 0 or 3 or when an analogous 

G cn op 

situation exists for C - 1 or 2. 

op 

Wagon-Wheel Cross-Sectional Geometry 

AT.c. If G = 2 or 3 and S = 3 Go to step Alia 

op op 

A17.a. Compute the radius of the fillet centers for a wagon-wheel grain. 



"b . Compute the angular location of the fillet centers 

0 = arc sin [(h + f)/A ] 

fw w fw 

A. 18. Compute the grain perimeters S„, and S,„ and initial cross-secitonal 

G1 (id 

areas and A^g^ using l = and a - and then 1=1^ and a - 

in the following equations . 
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Al8.a. If y- + f < a^ sin 6^ and y tan (a/2) <_ a, 

S G = n p la - 2y tan (a/2) + sin 6^ - f)/sin a - y cot a 

+ (f+y)(ir/2 + 6^) + + f + y)(ir/n p - 0^)3 

Go to step l8c 

If y + f > a^ sin 6^, y tan (a/2) _< a and 

if y < t 
J — wv 

Z 

S G * n p Uy+f)Wn p + arc sin (^j- sin 0^)3 

+ a. (ir/n - 0. )} Go to step l8c 

fw p fw 

if y > T w £ 

s g " n P (y+f){ v + arc sin L fzr sin V 

R 2 - - (y+f ) 2 

- arc cos I ] ) for S_ > 0 

2V y+f) 

- 0 for S_ < 0 Go to step l8c 

G G — 

Al8.b. If y tan (a/2) > a compute the coordinates, and y^' , 

(origin at fillet center) of intersections of burning slant sides 
with fillet arcs ; 

If a > 7t/2 , 

x^' = - a + y and 

y x ' = - [(f+y) 2 » x 1 ’ 2 ] Js 
If a < tt/2 

Q = - f + a tan a - y/cos a 

x^' = - Q tan a - I-Q 2 + (f+y) 2 sec 2 a] 2 /sec 2 a 
y^ 1 = x^' tan a + Q 



and the coordinates of the tips of grain points , 


y ' = - t, sin 0. and 
''o fw fv 

x 1 = - SL + y for a > tt/ 2 or 
o — 

x ' = (y ' - Q)/tan a for a < ir/2 

Then compute burning perimeters : 

If y + f <_ sin S^) 2 + x^' 2 ]** and 

s o = VI (x i' -* 0 ’> 2 + <V -y 0 ') 2 l h 

X T 

+ (f+y)l>/2 + - arc sin (^r)J 

+ (£„ + f+y)(iv/n - 0„ )} Go to step l8c 

IV p IV 

° r if y i T w> 

S r = n {[Cau* - x ') 2 + (y 1 1 - y ') 2 ]^ 

G p 1 o X O 

V 

+ (f+y) [tt/2 - arc sin ] 

R 2 - - (y+f) 2 

- arc cos ( )} for S_ > o 

s^iy+t) 0 

= 0 for S„ < 0 Go to step l8c 

G G — 

If y+f > IU^ sin 6^) 2 + aiia - 

if y < t , 

J — vw* 

Z 

S Q = n p {(y+f)Wn p + arc sin sin 0^)] 



-lU 


or if y > t , 
w 


fw 


S G = n p (y+f) {6 fw + arc Sin Sin 6 fw 


- arc cos [* 


R ~ 2 “ V ~ (y+f)2 


2 *fw (y+f) 


1> for S Q > 0 


S Q = 0 for S G £ 0 


Al8.c. If y £ 0, compute the (initial) cross-sectional area of the grain 
a Gs ■ Cl/2H»B 0 2 - sin 6^1 cos 0^ - (sin e^Koot o) 

- 2 U+f tan (a/2))/*,.^] - (ir-S^n^) J^ 2 - 2n^f[Jl 

+ ^ sin 0^) /sin a + - 6^) + U/n p 

+ ir/2 - esc a) f/2]} 

M9 * a - S G = S G1 + S G2 


A19 ‘ b - A Gs = A G1 + A G2 


Go to step A2k 


Initial Chamber Volumes 

A37.a. If y > 0 Go to step 50 

A37.b. If G =1, set L n . = 0 and A n =0 
op Gsi Gs 

37. c. If G =2, set L_ . - 0 and L_ . = 0 
op Gcx Gni 

37. d. Compute the approximate initial gaseous volume of the chamber. 


V c . - 1.10 {itD. 2 (L Qni ♦ L^.JA 


+ <« d 0 2 A - a Gs )l G s1 + Vrp "V A) ♦ v cir 


Go to step 50 


Note: The volume is approximate because an approximate relation is used for the 

volume of the thrust chamber passages and because an arbitrary ten percent 
additional allowance is made to account for slot volume, end volume and the 
effect of the convergent portion of the nozzle on the chamber pressure change 
characteristics. Any error thus introduced in the final volume and the 
volumetric loading density will be roughly one order of magnitude less than 
the error in initial volume. 
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Ignition Transients 


This analysis is treated in the computer program as a major sub- 
routine. The method used follows the general technique described by 
Sforzini and Fellows (Reference 2) for an igniter which is itself a small 
rocket motor. Features of the analysis include variation of the characteristic 
velocity of the main motor with chamber pressure. Chamber volume and surface 
area of the propellant grain are assumed constant during ignition, but the 
burning surface is established as a function of time in accord with the 
passage of the flame front which is assumed to spread at constant speed. Four 
changes are made in the referenced approach: 


1. A simple piecewise linear approximation to the pressure-time 
trace of the igniter (assumed to be a small SRM) is used 
rather than the polynominal approximation used in the 
reference. Comparison shows there is very little difference 
between computations by the two methods , and the program 
input is greatly simplified. 

2. The referenced analysis was modified to permit effects of erosive 
burning to be taken into account. This was accomplished by 
replacing the burning rate term Cpc 11 in equation 9 of Reference 

2 with the average burning rate as calculated by the methods used 
in the basic program of Reference 1. 

3. The igniter throat cross-sectional area is eliminated 

as an input variable by making use of an empirical criterion 
taken from Reference 3 for establishing the discharge mass 
flow rate of the igniter. 

k. Provision is made for including, as a program option, calcula- 
tions accounting for the effects of a main motor nozzle blowout 
plug. 

BO . Set t . = 0 

i g 

Bl.a. Compute the average mass flow rate of the igniter over the first half 

of the igniter burning time, excluding the pressure buildup phase 

(Time prior to tj^) 

« 

m. , s. ~ 0 ,2A*/g 
ig(av) 

This is an empirical relationship based upon Reference 3. 

Bl.b. Calculate the required igniter throat cross-sectional area 

' A*. -2m. , \ C*. /[2 P . -R . (t TO -t T1 )/2] 

ig ig(av) ig' L mig ng 12 II 
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B2. Compute the chamber pressure of the igniter 

B2.a. If t. < t T1 , 
ig — II* 

P . = P . t/t T - 

cig mig II 

B2.b. If t. > t_, and. P . > P 

ig II cig c 

P = P . - R . (t - t T , } 

cig mig rig II 

B2.c. If t. > t Tn and P . < P 

ig II cig — c 

P cig “ P c 

B3. Compute the estimated instantaneous mass flow rate of the igniter 


BU. 


If P . > P 

cig c 

m. = P . A*. /C* 

ig cig lg ig 

If P . < P or t. > t xo /2 

cig — c ig 12 

m. =0 (See discussion after KRIG, Section II) 
ig 

Solve the following equations for the chamber pressure of the main 
motor P c using a Runge-Kutta solution of the fourth order 


Bit .a. 

If 

P ■ 
c 


If 

P : 
c - 



”D 

Bk .b. 

r h 

- a 

Blt.c. 

C* 

= K 

BU.d. 

r 

= a 


big* setffi D =0 
big’ Set 

= A * P c /(K a + *b V 
c (i + r 2 J 2 /2)] n 

Vc 


n 


ujp A*X/C*A ) 0,8 (u. t. ) 0 ' 2 exp (fir p A C*/P A*X ) 
ph p # nn r s i .a n T) rm c 


eb 


pn 


fs ig 


P Pn 
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B5 . Compute 

P = P [l+r 2 J 2 /2] 
n c 

B6. Set P = P 

on c 

B7. Calculate C„, C„ , F, F , I , I m , I m using equations 

l * wloiouc -p vac’ * vac’ sp vac’ T* T vac 

20 through 2k with t. = t. . 

J -^*o 

B8. Set T. = l At. v Repeat step B2 through B7 

a . ^-6 * 

for k - 1,2,3. . . 

unt11 P <=(k) - P c(k-1) “ P o(k) /1000 
B9 . Calculate the estimated igniter propellant weight 

W ig(tct) = ®^ig(avj^ 5(t 12 " 



IV. THE COMPUTER PROGRAM 


This section contains the instructions for preparation and arrangement 
of the data cards for the extended program. A complete listing of the 
program statements is given folloved by flow charts for the two new sub- 
routines ; the ignition and the plotting subroutine for the "CalComp 663 " 
digital plotter. 

Data Card Usage 


The data formats have been established to allow the operator to look 
at the card and know which variables are represented on it. The format 
should be followed to insure correct reading of the inputs. The various 
data cards are as follows: (See Table IV-1 for the complete format) 



I. 

Number of configurations 


II. 

Initial zero values 

III 

a&b . 

User options (2 cards) 


IV. 

Propellant data 

V 

a&b . 

Motor geometry (2 cards) 

VI 

a&b . 

Performance data (2 cards) 


VII. 

Input, grain, etc. 

VIII 

a&b . 

Initial tabular inputs (2 cards) 

IX 

a&b . 

Tubular grain input (2 cards) 


X. 

General star data 


XI. 

Wagon -wheel data 


XII. 

Truncated star data 

XIII. 

Standard star data 


XIV. 

Termination port data 

XV 

a&b . 

Data for ignition transient calculations (2 cards) 

XVI 

a-e . 

Data for weight calculations (5 cards) 

XVII 

a&b . 

Tabular inputs (2 cards) 


NOTE: In the original report, tabular inputs and initial tabular 

inputs appeared under the same heading. 

Cards II-VII must accongpany each and every configuration even if only one 
parameter changes. If INPUT = 1 or 3, a number of sets of data cards 
XVII a&b are used. Consider a test run on four different configurations: 

a. A combination grain with only tabular inputs (no ignition 
or termination ports, but with weight calculations) . 
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b. A combination circular perforated and standard star grain with 
only equation inputs (No termination ports, no ignition, but with 
weight calculations). 

c. A truncated star grain with two termination ports using both 
tabular and equation inputs (No ignition * but with weight calcula- 
tions ) . 

d. A combination circular perforated and wagon-wheel grain using both 
tabular and equation inputs (No weight calculations or termination 
ports , but with ignition calculations ) . 

The correct order for the data cards is : 


Configuration a 


Configuration b 


1. 

I 


13 . 

11 

2. 

II 


i4. 

III a&b 

3. 

III a&b 


15 . 

IV 

4. 

IV 


16 . 

V a&b 

5. 

V a&b 


17 . 

VI a&b 

6. 

VI a&b 


18 . 

VII 

7. 

VII 


19 . 

IX a&b 

8. 

VIII a&b 

s 

11 

o 

20 . 

X 

9. 

XVI a-e 


21 . 

XIII 

10. 

XVII a&b 

(y=y-, ) 

22 . 

XVI a-e 

11. 

XVII a&b 

(y=yi) 



12. 

XVII a&b 

( 7 = 73 ) 




Configuration c Configuration d 


23. 

II 


37. 

II 


24. 

III a&b 


38. 

III a&b 


25. 

IV 


39. 

IV 


26. 

V a&b 


40. 

V a&b 


27. 

VI a&b 


4i. 

VI a&b 


28. 

VII 


42. 

VII 


29. 

VIII a&b 

O 

11 

>> 

43. 

VIII a&b 

(y=o) 

30. 

X 


44. 

IX a&b 


31. 

XII 


45. 

X 


32. 

XIV 


46. 

XI 


33. 

XVI a-e 


47. 

XV a&b 


3k. 

XVII a&b 

(y=y-, ) 

48. 

XVII a&b 

(y=y-, ) 

35. 

XVII a&b 

(y=yi) 

49. 

XVII a&b 

(y=yo) 

36. 

XVII a&b 

(y-yj) 

50. 

XVII a&b 
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Note that if a particular data card does not apply, no card is used. For 
example, XI is omitted for Configuration b. and XII for Configuration d. 

Program 'Listing 

A complete program listing is presented in Table IV-2. 

Flowcharts 

. Flowcharts of the ignition and plotting subroutines for the extended 
program are presented in Figures IV-4 and IV-5. 



TABLE IV-1 


DATA CARD FORMATS 


NUMBER OF CONFIGURATIONS TO BE TESTED =- 


0 .00 .00 .00 .00 .00 .00 .00.00 .00.00 .00 .00 .00.00 .00 .00.00 .00 .00 .00 .00 .0 


I GO-- IWO— 


IPO=- NUMFLT ( JJ )=- 


RHO— . A--.—— N=- . ALPHA- — . - BETA- MU-+-. E+— CSTAR- 


L= -. TAU- — . — DE- 


THETA--. ALFAN-- 


! LTAP- XT-—.— Z0- — . ■ 


DELTA Y=-. ZOUT- .— DPOUT- .— ZETAF- TB= HB— - 


RADER--. 


INPUT-— GRAIN-— STAR= 


NT- . ORDER-— COP=- 

I ! I 



TABLE IV-1 (Coat'd) 


ccl cclO cc20 cc30 ccUO cc50 cc60 cc70 cc80 


j \ i Villa i ! I 

Y= . — ABPORT=+-. E+ — ABSLOT=+- . E+ — ABN0Z=+-. E+ — 


lit! 

! ! VUIb ! i i 

APHEAD=+- . E+ — APHOZ=+- . E+ — VCIT=+-. E+ — 


i : I i 

! , ! •: iXa | I i 

j do= DI= . DELDI= . S= . THETAG=- . 

| | | ; j { i 

I I IXb i | 

i LGCI= LGNI= THETCK=-. THETCH=- . 

1 [ 


■ I I A 1 I 

NS= . LGSI= HP= . RC= . FILLET= . NN- 

} i , I 


1 i : ! Ai I ' j 

TAUWW= — . — Ll= — . — L2= — . — ALPHA1=-. ALPHA2=- . HW=— . 

■ ! : j | I ; ; 

i j ! XII ! i ! 

j RP= TAUS= . i ! i | j 


i i xm ; j 

THETAF=- . THETAP=- . TAUWS= i 

. , ! 


: I xiv 

LTP= . — DTP= — . — THETTP=-. TAUEFF=- . 


i ! . | avs . J t 

KA= KB=-. UFS= CSIG= PMIG= TI1=-.— 

i ■ , i i 


- i I 1 AV i 

TI2= — . — RRIG= DELTIG=- . PBIG= . 

i ' i t I 
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TABLE IV-1 (Cont'd) 


ccl 

? 

i 

PIPK: 

i 

; K2-- 


I LCC=- 


PSIA : 


KEH=. 


Y= 


cclO 


cc20 


DTEMP=- 


SYCNOM=- 


cc30 


cc^O 


cc 50 


cc 60 


[ I 

1 1 

' 

• 1 

1 

XVIa 



NSEG=- 


Kl=-.- 


KEN=-.— 


ABPOBT=+- . 


i 

APHEAD=+- . 


SIGMAP=- 


XVIb 


DCC=- 


SIGMAS= 


PSIC= 


HCN=- 


K2=- 

I 


XVI c 
SYMOM=- 


XVld 

PSIINS=- 


XVIe 


DLINER=- . • 


TAUL=- . - 


•E+- 


-E+- 


XVII a 

ABSLOT=+- . - 


XVI lb 
APNOZ=+-.- 


PSIS=- 


DELINS=- 


WA=- 


-E+- 


-E+- 


ABNOZ=+-.- 


ccTO 

1 


Nl=- 


DELC=- 


I 

-E+- 


cc80 

1 
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GE 


1 


TABLE IV-2 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


p 


c 

c. 

c 


c 

c 

c 

c 

c 


► 

c 

c 


**************** ********************* * ******* * ******* ******** *** ***** 


* SRM OcSIGM AND PERFORMANCE ANALYSIS * 

* PREPARED AT AUBURN UNIVERSITY * 

* UNDER MOD. NO. 8 TO COOPERATIVE AGREEMENT WITH * 

* NASA MARSHALL SPACE FLIGHT CENTER * 

* * 


* ANALYSIS BY R.H. SFORZINI * 

* PROGRAMMING BY J.M. LYON AND J.E. MURPH * 

* AEROSPACE ENGINEERING DEPARTMENT * 

* APRIL L973 * 

***** ***** ********** ** * * * ** ** * * & * * * *** * * * * * * ** *** * *** ***** ** * **** *** * 


INTEGER GRAIN 

real MGENtMDI $,MNOZ ,MN 1 , JROCK , N , L , MEL , ME , I SP , I T OT , MU , MASS ♦ I SPV AC 

REAL Nl,N2,NSEG,Kl,K2,KEH,KEN,NS,LCC,LTAP 

REAL M2 , MCB AR , ISP2»ITVAC,KA»KB » LAMBDA 

COMMON/ CON ST 1/Z W,AE » AT, THE TA , ALF AN 

C0IMM0N/C0NST2/CAPGAM,ME ,BOT, ZETAF ,TB,H8,GAM 

COMMUN/CONST3/S,NS,GRA IN 

COMMON /VAR I Al/Y,T,DELY,DELTAT,PONOZ,PHEAD,RNOZ, RH EAD, SUMAB , PHM AX 
COMMON/VARIA2/ABPORT.ABSLOT, ABNOZ,APHEAD,APNOZ ,OADY, ABP2 , A8N2 » ABS2 
COMMON/ VAR I A 3./ 1 TOT, I T VAC , JROC K , I SP , I SPVAC , MOI S , MNOZ , SG , SUMMT 
COMMON/ VAR I A4 /RNT, RHT, SUM2 ,Ri ,R2 , R3 ,RH AV E , KNAVE , RB AR , YB , KOUNT , T L 
COMMON/ VAR I A5/ABMAIN, ABTO , SUMDY , VC I , VC 
CO M.MON / I GN 1 /K A , KB , UF S » RHO , A ,N , L , PMI G , TI L ,T I 2 , C S I G 
COMMON/ IGN 2 /ALPHA, BETA , P9 I G, RR I G , DE LT I G , X , TOP, ZAP 
C0MM0N/PLUTT/NUMPLT(16) , IPO,NOUM,IPT,IOP 
DATA PI, G/3. 14159,32. 1725/ 

READ! 5, 500 I N RUN S 

****************************************************** *************** 


* READ IN THE NUMBER OF CONFIGURATIONS TO BE TESTED * 

********** *** *** i ***************************************************** 


500 FORMAT! 42X, 12 > 

I0P = 0 

DO 901 l -= 1 , NR UN S 
WRITF(6, 602) I 

602 FORMAT! lHl,42X, •CONFIGURATION NUMBER ',12) 

READ! 5,499) SUMDY , AN S , ZW , Y , T , DEL T AT ,RNGZ , RH EAD , SUMAB , PHM AX , SUM2 , IT 
IDT ,RHT,RNT,R1 , R2 , R3 , RH A VE , RNA VE , R8 AR , I T VAC , SUMMT 
**************** ************************** **** S'*** ******************* 


* SET INITIAL VALUES OF SELECTED VARIABLES EQUAL TO ZERO * 

* ***NOTE*** THESE VALUES MUST BE ZEROED AT THE BEGINNING OF * 

* EACH CONFIGURATION RUN * 

******** ***** *** ***************************************************** 


499 FORMAT! 22F3.1 > 

READ! 5,491) IGC,IWO 

READ! 5, 493) I PC » ( NUMPL TI J J ) ,JJ=1,16) 
********************************************************************* 

* READ IN THE USER'S OPTIONS * 


-2h~ 


>AGE 


2 


TABLE IV-2. (Cont'd) 


C 

c 

c 

c 

c 

c 

c 

►c 

►c 

►c 

>c 

>c 

►c 

>c 

l>C 

>c 

► 

►c 

►c 

►c 

►c 

►c 

>c 

>c 

>c 

>c 


►c 

►c 

>c 

>c 

►c 

>c 

►c 

>c 

c 

► 

► 

> 


c 

c 

c 

c 

c 

c 


* 


* 


* 

* 

* 

* 

* 

* 

A 

* 

* 

* 

* 

* 

* 


values for igo are 

0 FOR NO IGNITION TRANSIENT CALCULATIONS 

1 FOR IGNITION TRANSIENT CALCULATIONS 
VALUES FOR I WO ARE 

0 FOR NO INERT WEIGHT CALCULATIONS 

1 FOR INERT WEIGHT CALCULATIONS 
VALUES FOR IPO ARE 

0 FOR NO PLOTS 

1 FOR PLOTS OF EQUILIBRIUM BURNING ONLY 

2 FOR PLOTS OF IGNITION TRANSIENT ONLY 

3 FOR PLOTS OF BOTH IGNITION TRANSIENT AND 

EQUILIBRIUM BURNING 

VALUES FOR NUMPLTUJ) ARE (NOT REQUIRED FOR IPO-O) 

0 IF SPECIFIC PLOT IS NOT DESIRED 

1 IF SPECIFIC PLOT IS DESIRED 

CONTINUE 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


♦ 

* 

* 

❖ 


ORDER OF SPECIFICATION OF NUMPLTUJ) IS 

1 PHEAD VS TIME 

2 PONOZ VS TIME 

3 PHEAD AND PONQZ VS TIME 

4 RHEAD VS TIME 

5 R NO Z VS TIME 

6 RHEAD AND RNOZ VS TIME 

7 SUMAB VS TIME 

8 SG VS TIME 

9 SUMAB AND SG VS TIME 

10 F VS TIME 

11 F VAC VS TIME 

12 F AND F VAC VS TIME 

13 VC VS TIME 

14 SUMAB VS Y 

15 SG VS Y 

16 SUMAB AND SG VS Y 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


##*#*******###** X;***#*##********##***#*:*#^#*## **************** ******* 


491 FORMAT' 4X« lit 9X V 11) 

49 3 FORMAT { AX, U t 15X 1 16( 1 1 9 LX )) 

WR I TE ( 6 1 492 ) IGOtlWO 

WP IT E < 6 1 49 A ) IPO* ( NUMPLT ( J J) , JJ^l ,16) 

492 FORMAT* //, 2GX OPTIONS 1 ,/ t L3X f «IGQ* 1 , 1 1 , / , 13X , VI U 0* f ,Il) 

494 FORMAT* 13X, f I PO- 1 , 1 1 , / , 1 3 X , 1 NUMP LT ( J J ) * 1 , I 1 , 1 5( 1 , • , 12 > > 

READ( 5, 501 ) RHO* A, N, ALPHA f BETA, MU ,C STAR 

***$****$**$£*$£ ^^^ c ************************** ** * * ******************* 


* READ IN BASIC PROPELLANT CHARACTERISTICS * 

* $ 

* RHO IS THE DENSITY OF THE PROPELLANT IN SLUGS/IN*#3 * 

* A IS THE BURNING RATE COEFFICIENT * 

* N IS THE BURNING RATE EXPONENT * 
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GE 


3 


TABLE IV-2. (Cont'd) 


C 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


* ALPHA ANC BETA ARE THE CONSTANTS IN THE EROSIVE BURNING * 

* RELATION OF ROBILLARD AND LENIGR * 

* HU IS THE VISCOSITY OF THE PROPELLANT GASES * 

* C S TAR IS THE CHARACTERISTIC EXHAUST VELOCITY IN FT/SEC * 


i **************** ***** *** **** ****************** ******* ** ******** ****** 

501 FORMAT ( 4X,F8.6,3X,F6.4,3X,F5. 3, 7X , FA- 1 »6X , F5. I ,4X, E I I . 4 t 7X , F5. 0 ) 

WR ITE(6,603) RHO , A ,N, ALPHA , BETA ,MU,CSTAR 
603- FORMAT! //, 20 X ,* PROPEL LANT CHARACTERISTICS* ,/»13X»* RHQ=- *,F8.6»/,1 
1 3X , * A= * » F6 »4 , / , 13X * * N- • ,F5. 3 »/ 1 13X, • ALPHA* * , F4. 1 , / , 13X » ' 0ETA= ' 
2.F5. 1,/,13X,'MU* * »lPEL1.4t/«13Xt 'CSTAR= ',1PE11.4I 
READ! 5» 502 ) L ,TAU»DE*OTI » THETA .ALFAN.LTAP , XT , ZO 

***** ****** ***************************************************** 


* READ IN BASIC MOTOR DIMENSIONS * 

* * 

* L IS ThE TOTAL LENGTH OF THE GRAIN IN INCHES * 

* TAU IS THE AVERAGE WEB THICKNESS OF THE CONTROLLING GRAIN * 

* LENGTH IN INCHES * 

* DE IS THE DIAMETER OF. THE NOZZLE EXIT IN INCHES ■* 

* DTI IS THE INITIAL DIAMETER OF THE NOZZLE THROAT IN INCHES * 

* THETA IS THE CANT ANGLE OF THE NOZZLE WITH RESPECT TO THE * 

* MOTOR AXIS IN RADIANS * 

* AIFAN IS THE EXIT HALF ANGLE OF THE NOZZLE IN RADIANS * 

* LTAP IS THE LENGTH OF THE GRAIN AT THE NOZZLE ENO HAVING * 

* ADDITIONAL TAPER NOT REPRESENTED BY ZO IN INCHES * 

* XT IS THE DIFFERENCE IN WEB THICKNESS ASSOCIATED WITH LTAP * 

* ZO IS THE INITIAL DIFFERENCE BETWEEN WEB THICKNESSES AT THE * 

* HEAD AND AFT ENDS OF THE CONTROLLING GRAIN LENGTH * 


************** ** ******************* *** ******************************* 
502 FORMAT! 2X.F5. 0 , 5X, F 5 . 2 ,4X , F6. 2 , 5X , F6. 2 ,7X , F6.4 , 7X , F6 .4 , / , 10 X, 

I F6.2, 4X,F5. 2.4X.F5.2 I 

WRITE (6.60 4) L » TAU, DE , DTI » THE TA , ALF AN , LTAP , XT , ZO 
604 FORMAT ( / /, 20X , 'BASIC MOTOR DI MENSIONS* , / , 13X, • L= • , F5 .0 , / , 13X, 1 TAU 
1= ?, F5.2, /» 13X, »DE= '» 

21PE11.W, 13X , 'DTI-* , 1PE11.4,/,13X,'THETA* • , 1PE1 1 - 4 , / » 13X , • AL FAN= 
3 ' , 1PE11.4,/, 13X, ‘LTAP* •, 1 PE 1 1 . 4 ./ ,13X XT= » » I PE 1 1 .4, / , 13X , * Z0= 
4' , IP E 1 1 « 4 ) 

READ! 5 T 503 ) D ELT A Y , XOUT , DP OUT , ZE TAF ,TB ,HB »G AM , R ADE R 
**************** ***************************************************** 


* READ IN BASIC PERFORMANCE CONSTANTS * 

* * 

* OELTAY IS THE DESIRED BURN INCREMENT DURING TAILOFF IN INCHES * 

* XOUT IS THE 0 1 STANCE BURNED IN INCHES AT WHICH THE PROPELLANT 

* BREAKS UP 

* DPOUT IS THE DEPRESSURIZATION RATE IN LB/IN**3 AT WHICH THE * 

* PROPELLANT IS EXTINGUISHED * 

* ZETAF IS THE THRUST LOSS COEFFICIENT * 

* TB IS THE ESTIMATED BURN TIME IN SECONDS * 

* HB IS THE ESTIMATED BURNOUT ALTITUDE IN FEET * 
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* * 


AGE 


TABLE IV-2. (Cont'd) 


C * GAM IS THE RATIO OF SPECIFIC HEATS FOR THE PROPELLANT GASES * 

C * RADER IS THE RAOUL EROSION RATE OF THE NOZZLE THROAT IN * 

C * INCHES/SEC * 

£ #***###**$***$*# ****** ******* ************* ************ ******* ******** 
503 FURMATI 7X , F5. 3 ♦ 6X, F 7. 2 , 7X , F 7. 2 , 7 X ,F A. 3 ,4X, F5 . 1 , 4X , F7. 0 , 5X , F5 .3 , / , 8 
1X,F6.4) 

WRITE (6,606) DEL TA Y , XQUT , DPOU T , Z6 TAF , TB ,HB , GAM , RADER 
606 FORMAT! //, 20X »• BASIC PERFORMANCE CONSTANTS* , / » 1 3X , * DELT AY= SF5.3, 
1/, 13X, »XOUT= *,F7«2,/,13X,* DP OUT » SF7.2,/ ,13X , • ZETAF= *,F5.3,/,13 
2X,'TB= SF5.lt/tl3X, *HB* • ,F7.0, / , 13X, *GAM= * , F5.3 ,/ , 13X , * RADER= ' 
3, FEU) 

► N0IJM=0 

► I PT = 0 
MN 1= . 65 
ME 1=7.0 
Z = ZO 

S = 0.0 
N S = 0 . 0 
KGUN T = 0 
A8MAIN=0.0 
ABTO=0.0 
DELy=DELTAY 
T OP = GAM + 1 • 

B0T=GAM-1. 

7AP=T0P/< 2.*B0T) 

CAPGAM= SORT (GAM )*( 2. /TCP )**ZAP 
AE=PI*0E*0E/4. 

1 IF(XT.LE.O.O) TL=0 . 0 
IF(XT.LE.O.O) GO TO 40 
TL=< Y-TAU+XT+Z/2. )*LTAP/XT 
I F ( TL .LE.O.O) TL=0 . 0 

IF( TL .GE.LTAP ) TL=L TAP 
40 0T=DTU2 .*< RAO£R*T) 

AT=Pl*DT#DT/4. 

CALL AREAS 
IF(Y.Lt.O.C) VC=VCI 
IF(AOS( ZW) .GT.O.O) GO TO 20 
IFISUMAS. LE.O.O) 60 TO 31 
X= { ABP0RT+A8SL0T ) / SUMAB 

90 MNCZ=AT*X/APN0Z*(2.*(1.+80T/2.*MN1*MN1)/T0P)**ZAP 
IF( A9S( MN0Z-MN1) .LE. 0.002) GO TO 2 
MN l=MNQ Z 
GO TO 90 

2 VNUZ=GAM*CSTAR*MNOZ*SQRT ( ( ( 2. /TOP) ** ( TOP/ BOT ) ) / ( l . +80T/2 .*MNOZ *MNO 
1Z ) > 

PR AT = ( 1 .+BGT/2.*MNQZ*MN0Z) ** (-GAM/BGT) 

JR0CK=AT/APN0Z 

IF( IGO.EQ.G.OR. Y. GT.O.O) GO TO 900 
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GF 


5 


TABLE IV-2. (Cant'd) 


REAUt 5,57) KA , KB , UF S,C SI G ♦ PMI G ,T 1 1 ,T 1 2 , RR1 G , DE L T IG , P8IG 
C ******** **************************************************** ********* 
C * READ IN VALUES REQUIRED FOR IGNITION CALCULATIONS * 

C * ***NOTF*** NOT REQUIRED IF I G0=0 * 

C * 

C * KA AND K6 DEFINE THE CHARACTERISTIC VELOCITY IN FT/SEC 

C * CSTA * KA «• K8 * PRESSURE 

C * UF S IS THE FLAME-SPREADING SPEED IN IN/SEC 

C * C S IG IS THE CHARACTERISTIC VELOCITY OF THE IGNITER IN FT/SEC 

C * PM I G IS THE MAXIMUM IGNITER PRESSURE IN L8S/IN**2 


C * Til IS THE TIME OF MAXIMUM IGNITER PRESSURE IN SECONDS * 
C * T 1 2 IS THE TIME! IN SECONDS) FOR THE IGNITER PRESSURE TO * 
C * DROP TO 10 PER CENT OF MAXIMUM VALUE(PMIG) * 
C * RR IG IS THE AVERAGE REGRESSION RATE OF THE FIRST HALF OF THE * 
C * IGNITER PRESSURE TIME TRACE IN LBS/ I N**2/S EC * 
C * DELTIG IS THE TIME INCREMENT FOR IGNITION TRANSIENT * 
C * CALCULATIONS IN SECONDS * 
C * PBIG IS THE BLOWOUT PRESSURE OF THE MAIN MOTOR BLOWOUT PLUG * 
C * IN LBS/IN**2 * 


C **************** * ******** ********** **#*#**#*#***#**$#* ******* ******** 
97 FORMAT ( 3X , F7. 1 , 5X , F5. 3 , 6X , F 7. 1 ,7X ,F 7. 1 , 7X , F7. I ,6X , F5. 3 , / ,4X » F5.Z, 

I 7X,F7. L,5X, F5.3, 7X,F7.3) 

WR IT E( Of 8.42) KA, KB, UF S, C SI G, PMIG , Til, TI2,RPIG, DELTIG, PBIG 
642 FORMAT ( 20X , * I GN I T ION CONSTANTS* ,/ ,13X, » KA= » , F7 . 1 , / , 1 3X , * K8= *, 

1 F5 • 3, / » 13X, ' UF S- ',F7.l,/,13X,'CSIG= • ,F7. 1 ,/ ,13X PMIG= *, 

1 F7.1,/,13X, *TIl= » ,F5.3,/fl3X, *TI2= * ,F5 . 2 . / , 13X , ■ RR I G= *, 

1 F7.lt/, 13X, *DELT1G= • ,F 5. 3 , / , 1 3 X , ■ PBIG= *,F7.3,//) 

900 IF( IWG.FQ.C.OR.Y.GT.O.O) GO TO 832 

READ! 5,600) P I PK , DTE MP , SI GMAP , SI GMAS, N1 , N2 , S YCNOM , DCC, PS I C, DELC, LC 
1C,NS£G,HCN, SYNNCJM,PSIS,P$IA,K1,K2,PSI INS ,DEL INS , KEH , KEN, DL INER, TAU 
2L,WA 

C **************** ***************************************************** 


C * -READ IN BASIC PROPERTIES REQUIRED FOR WEIGHT CALCULATIONS * 
C * ***NOTE*** mot REQUIRED IF I WQ=0 * 
C * * 
C * P I PK IS THE TEMPERATURE SENSITIVITY COEFFICIENT OF PRESSURE * 
C * AT CONSTANT K * 
C * DTEMP IS THE MAX EXPECTED INCREASE IN TEMPERATURE ABOVE * 
C * CONDITIONS UNDER WHICH MAIN TRACE WAS CALCULATED IN * 
C * DECREES FAHRENHEIT * 
C * SIGMAP IS THE VARIATION IN PHMAX * 
C * SIGMAS IS THE VARIATION IN CASE MATERIAL YIELD STRENGTH * 
C * N1 IS THE NUMBER OF STANDARD DEVIATIONS IN PHMAX TO BE USED * 
C * AS A BASIS FOR DESIGN * 
C * N2 IS THE NUMBER OF STANDARD DEVIATIONS IN SY TO BE USEO AS * 
C * A BASIS FOR DESIGN * 
C * SYCNOM IS THE NOMINAL YIELO STRENGTH OF THE CASE MATERIAL * 
C. * IN LBS/INCH * 
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DCC IS THE ESTIMATED MEAN DIAMETER OF THE CASE IN INCHES * 

PSIC IS THE SAFETY FACTOR ON THE CASE THICKNESS * 

DELC IS THE SPECIFIC WEIGHT OF THE CASE MATERIAL IN LBS/IN**3 * 
LCC IS THE LENGTH OF THE CYLINDRICAL PORTION OF THE CASE * 

INCLUDING FORWARD AND AFT SEGMENTS IN INCHES * 

NSEG IS THE NUMBER OF CASE SEGMENTS * 

HCN IS THE AXIAL LENGTH OF THF NOZZLE CLOSURF IN INCHES * 

SYNNCM IS THF NOMINAL YIELD STRENGTH OF THE NOZZLE MATERIAL * 
IN LBS/INCH * 

CUNT INUE 

PSIS IS THE SAFETY FACTOR ON THE NOZZLE STRUCTURAL MATERIAL * 

°SIA IS THE SAFETY FACTOR ON THE NOZZLE ABLATIVE MATERIAL * 

K1 AND K 2 ARE EMPIRICAL CONSTANTS IN THE NOZ7LE WT. EQUATION * 
PSIINS IS THE SAFETY FACTOR ON NOZZLE INSULATION * 

DELINS IS THE SPECIFIC WEIGHT OF THE INSULATION IN LBS/IN**3 * 
KEH IS THE EROSION RATE OF INSULATION TAKEN CONSTANT * 

EVERYWHERE EXCEPT AT THE NOZZLE CLOSURE IN IN/SEC # 

KEN IS THE EROSION RATE OF INSULATION AT THE NOZZLE CLOSURE * 
IN IN /SEC * 

DLINER IS THE SPECIFIC WEIGHT OF THE LINER IN LBS/IN**3 * 

TAUL IS THE THICKNESS OF THE LINER IN INCHES * 

WA IS ANY ADDITIONAL WEIGHT NOT CONSIDERED ELSEWHERE IN LBS * 

*#**V*****$>Sl##ik# **$$**#* $*4^ ^ 

600 FORMAT! 7X » F 5* 4. 9X, F6.2, iOX *F5*3, 10X *F5* 3 ,6X * F5. 2»/ ,5X* F5 . 2 , 10X, FIO 
1.2f7X»P6.2»8X»F5*2»8XtF5.3»/f6X*F7.2»8X»F3.0»7XiF4.l»10X»FI0.2»8X, 
2F5.2t/»/X»F5. 2»6X«F6«4»6X»F6*4,10X»F5.2*10X,F6.4,/*6XpF6.4t7XfF6.4 
3,10X,F6.4,£X,Fc.4«6X,F7.2) 

WR ITE I 6, 610 ) PIPK T OTEMPp SIGMA P»SIGMAS»N1 »N2 t S YC NOM» DCC * PSICpDELC»L 
ItC tNSEG»HCN» S YNNOM *PSIStPSIA,KL»K2 f PSIINStDELINS,KEH»KEN,DLINFR,TA 
20L » WA 

610 FORMAT! 20X, 'INERT WEIGHT INPUT S , t/»13X»*PIPK= * « IPEi 1.4, / , 13X, * 

1 DT EM P= '»1PE11.4»/ » 13X» * SIGMAP- * *1 PE 11* 4*/ *13X»* SIGMA S= ‘.1PEL1.4 
2»/» L3X, »Nl= ' * 1PE11.4,/, 13X,»N2= ' » 1PE11.4 , / , 13X , * SYCNGM= • *lPEli. 
34»/» 13X, 'DCC= ' » lPE11.4t/«13X, *PSIC= • , 1 PE 1 1 . 4 , / , 1 3X , 1 DELC - ' t lPEl 

41.4. /,13X,«LCC= lPE11.4,/ f 13X t *NSEG= % l PE 1 1 . 4 ,/ , 13X ,» HCN= » , 1PE 

511 .4, /» 13X, 'S YNNOM = » , lPEll.4f / ,1 3X » • PS I S= * , 1 PEI 1 ,4» / , I3X * • PS I A= 
6* , lPEll.4f/»13Xt 'KI* * tlPEl 1.4 »/, 13 X, 'K2= • ,1 PE11 . 4 ,/ r 13X , • PS 1 1 NS = 
7 * » 1PE 11. 4* /, 13X » ' DEL INS= 1 * 1PE1 1 . 4 ♦/ , 1 3X » * KEH= • » 1PE1 l .4» / » 13X » 1 K 
8EN= ' • 1PE11.4*/ f 13 X » 'DLINER® • , 1 PEI 1. 4 , / , 13X , • T AUL= ' , 1 PE 1 1 .4 f / . 13 
9X, *WA» S iPEl 1.4) 

832 SUMY A =D ELY* (A BP2+ABN2+A8S2 ) 
lE(Y.EQ.O.O) SUMYA = 0.0 
VC=VC+SUMYA 
IF(Y.GT.G.C) GO TO 11 

PONOZ*! A*RHG*CSTAR*$UMAB/AT)**( 1. /{ l.-N} ) * ( 1.+ ( CAPGAM*JROCK J **2/2 . 

1 )**< N/( 1.-N ) ) 

PON=PONOZ 

MDIS=AT*PONOZ/CSTAR 


C * 
C * 
C * 
C * 
C * 
C * 
C * 

c * 

C * 

c * 
c * 
c * 
c * 
c * 
c * 
c * 
c * 
c * 
c * 
c * 
c * 
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P2=P0N0Z 
PQNOZ 2= PONGZ 
PN0Z=PRAT*PONOZ 

P4=2 . *MDIS*VNOZ/(APHEAD+APNQZ )+PNQZ 
I F ( GRAIN .60 .3 ) P 4=MD I S*VNOZ/APNOZ+PNOZ 
5 PNOZ=PRAT*PONOZ 

PHEAU=2 .*MD I S * VNGZ / l APHE AD+AF NQZ ) + PNOZ 

I F < GR A I N • EQ .3 ) PHE AD=MOI S*VNOZ/APN0Z+ PNOZ 

kHEAO=A*PHEAD**N 

Z IT=MDI S#X/APNOZ 

RN 1=R HEAD 

PHEADZ=PHEAD 

3 RNOZ =RN l-( ( RN 1-A*PN0Z**N-ALPHA*ZIT#*. 8Z (L#*.2*EXP< SET A*RNi *RHQ/ Z IT 
I) } )/( l.+AtPHA*ZIT**.8*BeTA*RH0/ZiT/a**.2*EXP(8ETA*RNl*RH0/Zm ) )) 

IFtAbSIRNl-RNOZ) .LE. 0.002) GO TO 4 
RN 1=RN0Z 
GO TO 3 

4 A V E 1 = ( R HE AD+R NO Z ) /2 • 

IFtY.GT .0.0) GO TO 7 

RN2=RN0Z 

RH2=RFEAD 

PON J=PONGZ 
OP CDY = 0.0 
AV 62 = AV E 1 

7 RNAVE=( RN0Z+RN2) /2. 

RHAVE=( RHEAD+RH2 )/2. 

MGEN=RH0/2.*( ( RNOZ +R HEAD ) * ( ABP0RT+A9SLGT) f2. *A*PQNQZ**N*ABNOZ ) 

DROY={ AVE I-AVE2 ) /OELY 

RBAR=(AVFl+AVE2)/2. 

GMAX=1.C02*MD IS 
GM IN=0.998*M0IS 
IF(Y.GT.O.C) GO TO 12 
GMAX=1.00L*MD IS 
GM IN=0.999*MD I S 

IF ( MGEN ,GE .GM IN.AND.MGEN.LE.GMAX) GO TO 6 
M01S*MGEN 

P0NOZ=MD I S*CS TAR/AT 
GO TO 5 

6 RE=2.*MDIS*X*L/< < APNO Z+APHE AD ) *MU) 

17 ME=SQRT(2./B0T*< T0P/2.*( AE*ME l/AT ) ** l 1 . /ZAP )-l . ) ) 

IP UBS( M6-ME1) .LE. 0.002) GO TO 9 
ME I=M6 
GO TO 17 
9 CONTINUE 

IF ( IGO.NE.O.AND.Y.LE.O.O) CALL IGNITN 
IF (Y .LE .0.0 ) WR l TE ( 6» 1 01 ) RE 

101 FO RM ATI / //i 33X» * if*#*####*****#**# ************* » / »33X»***** EQU 1 L I B 
1R I UM BURN ING ***** , /» 33X» **************** *************** »// f 30X t 
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l 'INITIAL REYNOLDS NUMBER* *tlPE11.4) 

POM J - PONU7. 

CALL OUTPUT 

10 IF(Y.LE..O 5*T AU ) GO TO 16 

SINK 1=VC/(CAPGAM*CSTAR)**2*RBAR*DPC0Y/12. 

MA SS = . 0 1*MD I S 
AN$4=Y+10.C*DELTAY 
IFIKOUNT.GT.G ) GO TO 16 

1F( ABSl SINKU ,LE. MASS. AND. ANS4.LE.ANS-XT) GO TO 18 
GO TO 16 

18 DELY=iO.#DELTAY 
GO TO 5 5 

16 OELY=DELTAY 

55 DELTAT=2.*0ELY/(RHAVE+RNAVE) 

Z = Z ♦DEL TAT* (RNAVE-RHA VE ) 

Y=Y+DEL Y 

T = T+OEL TAT 

SUM2=SUMAB 

RN2*RN0Z 

RU2=RFEA0 

AVE2=AVEi 

GO TO 1 

11 MOIS=AT*PONOZ/CSTAR 
GO TO 5 

12 DP COY - ( PHEA02+PCN0Z2) / (RNAVE+RHAVE) *DRDY+ ( PHEAD2+P0N0Z2 ) / I ( ABP2+A0 
1N2+ABS2 J *2 . )* DAOY 

IF ( ABSI DPCUY) .GE.DPQUT.OR.Y.GE.XOUT) GO TO 25 

SINK i = V C/ ( CAP GAM*CSTAR >**2*RBAR*DPCDY/12. + < PHEAD2 + PONOZ2 >/2.*< RNAV 
1E + RHAVE >/2.*{ ABP 2+ABN2 *-ABS2)/(12.* ( CSTAR*C APGAM ) **2 ) 

ST UFF=M GEN-SI NK 1 

IFISTUFF.GE.GMIN. AND. STUFF. LE.GMAXI GO TO 14 
MD IS* STUFF 
PUNOZ=MOI S*CSTAR/AT 
GO TO 5 
14 P1=P0N0Z 
PON J =PQNOZ 
PCN0Z2=(Pl+P2)/2. 

P2=P0N0Z 

P3=PHEAD 

PHEA02=IP3+P4)/2. 

P4=PHFAD 

AN S = T AU- AB SI Z /2 • ) 

IFIY.LT.ANSl CALL OUTPUT 
IF(Y.LT.ANS) GO TO 10 

19 ZW=7 
YW=Y 

SUMBA*SUMAB 

P1=P0N0Z 
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KH2=RHEAU 
RN2=RNGZ 
RAVE=AVE1 
A8MA I N = SUMAB 
At3 TQ=0. 0 
WRITE(6,5i > 

51 FORMAT! //, 37X ,'*************************»,/ t 37X,« **** TAIL OFF 8EG 
LINS ****», / t 3 7Xt ******************* **#****• ,/) 

20 ANS2 = TAU+ABS< ZW/2. ) 

KOUNT =KOUNT +1 
DELYW=DELTAY 

OY 2=DEL YW 
IF(ZW) 32.32,33 

32 IF(Y.LT.ANS2.AND.ABS(Zk).GT.0Y2) GO TO 21 L 
SUMA8=A8MA IN 

GO TO 31 

211 SUMDY = SUM QY+DELY W 

SUMAB = ( U+SUMDY/ZW-DEIYW/<2.*ZW) > *ABTO- ( SUMDY/ ZW-DELYW/ ( 2 .*ZW > ) *AB 

IMA IN 
GO TO 31 

33 IF(Y.LT.ANS2.AND.ZW.GT.OY2> GO TO 21 
SUMAB=ABTO 

GO TO 31 

21 SUMOY=SUMUY+OELYW 

SUMAB = ( l.-SUM0Y/ZW+0ELYW/(2.*ZW) )*A8MAIN+< SUMDY/ZW-DELYW/ ( 2 .*ZW ) >* 

1 Art TO 

31 IF( SUMAB. IF.O .0) PONOZ=PONQZ /2. 

IF(SUKAB.LE.O.O) GO TO 25 

PGNOZ = ( A*RHO*C STAR*SUMAB/AT ) *# ( l * / < l.-N ) ) 

M D I S = AT *P ONOZ /CSTAR 
ARAVE = < SUMAB+SUMBA) /2. 

SUMYA=DELY*ABAVE 
VOVC + SUNYA 

CA DY = { SUMAB-SUMBA}/DELY 
P3AR=(P 1 + PONOZ )/2. 

SUM«A=SUMAB 

22 OPCOY=PBAR/( 1 .-N ) * l „ /ABA VE*DAD Y 
IF{PQNOZ.LE.30.0) GO TQ 25 
■RNOZ = A*PCNOZ**N 

RHEAD=RNOZ 

R 3AR = (RHEA0+RAVE )/2. 

MGCN=RHO#( RNOZ+RHEAD) /2.*SUMAB 
GMAX=1.002*MDI S 
GMIN=0.598*NDIS 

S INK l=VC/( CAP GAM*C STAR) **2*RBAR*DPCDY/12 . 

STUFF=MGEN-SI NKi 

IF(STUFF.GE .GMIN. AND. STUFF. LE.GMAX1 GO TO 23 
M D I S = ST UFF 
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PUNOZ =PCN J +DP C0Y*Q£LY 

IF ( PONOZ .LE*0 .0 ) PONO Z=0. 0 

P6 AR= ( P 1+P0N0 7 ) /2. 

GO TO 22 

23 RHAVE={RF2+RHEAD)/2. 

RN AVE=( RN2+RNOZ ) /2. 

RH2=RFEAC 

RN2=PN0Z 

PHEAD=PONOZ 

KAVE*RHfcAD 

P 1-PflNQ Z 

Pl>NJ=PONOZ 

IF ( A B S ( DPCDY) .GE.DPOUT) GO TO 25 
I F ( Y .GE *XOUT) GO TO 25 
CALL OUTPUT 
GO TO 10 

25 SUMAB=0.0 
RHEAU=0.0 
RN OZ -R HE AD 
PHEAD=PONUZ 
WRirE<6»3I8> 

318 FOKMATl //, 33X , ****************** ********** * ****** , / ,33X, ' BEG I 
IN HALF SECOND TRACE ***# * ,/ f 33X» * ********************************* 
I* > 

CALL OUTPUT 
TIME=T 
DELTAT= .5 
T I M=T IM E + 5. 

PHT =PH£ AD 
PONT = PONOZ 
SG=0 . 0 

29 T= T + DEL TAT 

PHEAO=PHT/EXP { CAPGAM**2* AT*C STAR / VC* ( T-T I ME ) *12. ) 

PONOZ =P HE AD 
MO I S= PONOZ *AT /CSTAR 
Y=Y + « 5*RHEA0 
CALL OUTPUT 

26 I*=(T.LT. TIM. AND. PHEAD.GE. 30.0) GO TO 29 
100 IriP 1 = G* SUMMT 

wP2=RH0*( VC-VCI )*G 
wP=( WPI+WP2I/2- 
ISP= I TOT /WP 
1SPVAC*ITVAC/WP 

Wft ITE(6 t 102) WPl*WP2,WP f PHMAX , I SP , I SPVAC , I TOT ♦ I TV AC 
102 FORMAT! ///, 13X, • WP1* ' , IPE 1 L. 4, / , 13 X » ♦ WP2= • , 1 PEI 1,4,/ . 13X , ' WP = '* 
11PE11.4,/»13X, * PHMA X = • , 1PE11.4,/ ,13X ♦ M SP= ',1PE11.4 t/» 13X,»ISPVA 
2C- ' * 1PE11.4* /*13X, 'ITOT* ' » 1P61 1 . 4 , / , 1 3X, ' I TVAC* •tlPE11.4) 
LAMBDA=( VC-VC D/VC 
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WRITE(6,103) VCItVC, LAMBDA 

103 FORM AT ( 13X , 4 VC I* * , 1PE 1 1 . 4 , / , 1 3X , * VCF= 4 , 1 PE 1 1 .4, / , 13X , • LAMBDA^ 

1 1PE11.4) 

1F( IWO.EQ.C) GO TO 903 

PMEOP=PHMAX*( l.+Nl*$IGMAP )*£XP(PI PK*DTEMP) 

SY C=S YCNUM # ( l • -N2# S IGMAS ) 

TAUCC= p SIC*PM EGP*DCC / ( 2»* SYC ) 

WCC=Pi*TAUCC*DCC*DELC*LCC*(i.+(NSeG-l. ) * (40. *TAUCC/LCC ) ) 
TAUCD=TAUCC/2. 

WCH=2.5* p I/2.*DCC*DCC*TAUCD*DELC 
WCN=4.5*PI/2.*DCC*HCN*TAUCD*D£LC 
WC=WCC+V>CH+WCN 
EP S IL-AE/AT 
DT=2.*SQRT< AT/PI ) 

WN=K 1*DT*0T /( 1.+.5*SIN(ALFAN) )*( { E PSI L-SQRT ( EPS I U ) * PMEQP*DT*PS IS/ 
1SYNN0M+K2*T*PSIA) 

WINSaT*PSl 1NS*DELINS4DCC*PI*( KEH* ( DCC*. 40+ ( S+NS ) *T AU/2 • +0. 15/ 

IPS IINS*(LCC-TAU*l S+NS) ) ) +KEN*. 80*HCN) 

WL=TAUL*DL INER*PI*DCC*(DCC/2.+LCC+HCN) 
wM = wC +WN+W I NS +WL + W A + WP 
Z E TAM =WP /WM 
RATIO=ITUT/WM 
WR I TE ( 6 » 605 ) 

605 FORM AT ( ///, 2QX, 'MOTOR HEIGHT CALCULATIONS') 

WR I TE ( 6 » 6 01 ) PMEOP, TAUCC ,WC,WN,WI NS , WL » WM , ZET AM , RAT I Q 
601 FORMATl 13X, 'MAX EXPECTED PRESSURE- 4 , 1 PEI 1 . 4 , / , 13X , ' CYL I NOR I CAL CA 
1SE THICKNESS^ ' , IP E 1 1 . 4 , / , 1 3X , 4 C A SE WT= ', 1PE 1 L . 4 ,/, 13X NOZZLE WT 
2= ', 1PE11. 4,/, 13X, 'INSULATION WT= » » 1 PE 1 1 - 4 , / , 13X , 4 LI NE R WT = 4 , 1PE 
311.4,/, 13X, 'TOTAL MOTOR WT= 4 , 1PE 1 1 .4 , / , 13X , 4 Z£ TAM= 4 , 1PE11 .4, / , 1 3 
4X, 'RATIO OF ITOT TO WM= 4 ,IPE11.4) 

903 CONTINUE 
NDUM=1 

IF( IPO.NE.O.AND. IPO.NE *2 ) CALL OUTPUT 
901 CONTINUE 
► IF(IOP.NE.O) CALL PLQ T ( 0. 0 , 0. 0,999 ) 

STOP 

END 
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SUBROUTINE AREAS 

C ********************************************************** **** *#*# $** 
C * SUBROUTINE AREAS CALCULATES BURNING AREAS AND PORT AREAS FOR * 

C * CIRCULAR PERFORATED (C.P.) GRAINS AND STAR GRAINS OR FOR A * 

C * COMBINATION OF C.P. AND STAR GRAINS * 

C ***** #4 ***# ***###****#* #**♦***# *#**#***#>(<**#** #4* *<!*#* ###*#***** 

integer star.grain.oroer.cop 

REAL MGEN.MDI S. MNOZ , MN 1 , JROCK * N . L » ME l » ME . I S P t I T OT , MU » MASS » I SPVAC 

REAL LGCI.LGN I.NS.NN.NP.LGSI » NT , LTP » LGC . LS , LF 

REAL M2 » MDBAR. » I SP 2 * I T VAC , LI ♦ L2 . LF W . LF WSQD 

COMMON/CON ST1/ZW.AE. AT .THETA .ALFAN 

C0MM0N/CQNST3/S.NS. GRAIN 

COMMON/ VAR I A1 / Y » T » DEL Y .DELTAT , PONQ Z . PHEAD.RNOZ , RHE AO, SUMAB , PHMAX 
COMMON/ VAR IA2/ABPORT, ABSL0T y A8N0Z .APHEAD.APNOZ ,DADY, A8P2, ABN2, ABS2 
COMMON/ VAR I A3 / 1 TOT , IT VAC , JROC K , I SP y l SP VAC y MDI S y MNOZ y SG y SUMMT 
COMMON/ VAR I A4/RNT,RHT»SUM2»R1 , R2 » R3 ,RH AVE , RNAVE ,RB AR » Y B , KOUNT » T L 
COMMON/ VAR I A5 /ABMA IN » ABTO » SUMDY » VC I » VC 
DATA PI/3.14159/ 

A6 PC- O. G 
ABNC-O. 0 
A G SC = 0 . 0 
ABPS-=0. 0 
ABNS -0. 0 
AiSS=0.0 
DABT^O.O 
SG-0 . 0 
ANUM=PI /4. 

PI D2=P 1/2. 

RNT=RNT ♦■RNOZ* DELTAT 
RHT=RHT+RHEAD*OELTAT 
IF(Y.LE.O.O) AGS=0. 0 

1 K=0 

lFIABS(ZW) .GT .0.0) K=1 
Y8=Y 

IF(K.EO.i) Y= Y8- SUMO Y/2. 

2 IF(K.EQ.Z) Y=YB+AQS(ZW) /2.-SUMDY/2. 

I^IY.Lt .0.0) READ! 5 y 500) INPUT»GRAI N»STAR, NT , ORDER. COP 

C # # #!(£$*#***# ifir It:#########** its ## Xofe #*&***# #***::<<:***>)«*#::* *** 

C * READ THE TYPE OF INPUT FOR THE PROGRAM AND THE BASIC GRAIN * 

C * CONFIGURATION AND ARRANGEMENT * 

C * VALUES FOR INPUT ARE * 

C * 1 FOR ONLY TABULAR INPUT * 

C * 2 FOR ONLY EQUATION INPUTS (EQUATIONS ARE BUILT * 

C * ' INTO THE SUBROUTINE) * 

C * 3 FOR A COMBINATION OF 1 AND 2 * 

C * VALUES FOR GRAIN ARE * 

C * l FOR STRAIGHT C.P. GRAIN * 

C * 2 FOR STRAIGHT STAR GRAIN * 
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C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

L 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


* 

* 

* 

* 

* 

$ 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


3 FOR COMBINATION OF C.P. AND STAR GRAINS * 

VALUES FOR STAR ARE (WAGON WHEEL IS CONSIDERED A TYPE OF * 

STAR GRAIN IN THIS PROGRAM) * 

0 FOR STRAIGHT C.P. GRAIN * 

1 FOR STANDARD STAR * 

2 FOR TRUNCATED STAR * 

3 FOR WAGON WHEEL * 

VALUES FOR NT ARE * 

0 IF THERE ARE NO TERMINATION PORTS * 

X WHERE X IS THE NUMBER OF TERMINATION PORTS * 

VALUES OF ORDER ESTABLISH HOW A COMBINATION C.P* AND STAR * 

GRAIN IS ARRANGED * 


1 IF DESIGN IS STAR AT HEAD FND AND C.P. AT NOZZLE * 

2 IF DESIGN IS C.P. AT HEAD END AMD C.P. AT NOZZLE * 

3 IF DESIGN IS C.P. AT HEAD END AND STAR AT NOZZLE * 

4 IF DESIGN IS STAR AT HEAD END AND STAR AT NOZZLE * 

***NGTE*** IF GRAIN*lf VALUE OF ORDER MUST BE 2 * 

***NOTE*** IF GRA I N =2 * VALUE OF ORDER MUST BE 4 * 


CONTINUE 

* VALUES FOR COP ARE (APPLICABLE TO C.P. GRAINS ONLY) * 

* 0 IF BOTH ENDS ARE CONICAL OR FLAT * 

* 1 IF HEAD END IS CONICAL OR FLAT ANO AFT END IS * 

* HEMISPHERICAL * 

* 2 IF BOTH ENDS ARE HEMISPHERICAL * 

* 3 IF HEAD END IS HEMISPHERICAL AND AFT END IS 

* CONICAL OR FLAT 


«****«)!< *«#«#**#* #*«*****♦ * #$«*$*♦* * 


300 FORMAT! 9 X* I 2 , 9 X . I 2 t 8 X ♦ I 2 » 6 X , F 4 . 0 , 9 X , 1 2 , 7 X , 1 2 ) 

IF(Y.LE.O.G) WR I TE ( 6 » 6 C 7 ) 

607 FORMAT! //, 20 X , 'GRAIN CONFIGURATION') 

IF(Y.LE.O.O) WRITE! 6 » 600 ) INPUT* GRAIN. STAR. NT tCRDER.COP 
600 FORMAT! 13 X, ' INPUT* ’ » I 2 . / , 1 3 X , « GR A 1 N= * , I 2 , / . 1 3 X, • ST AR= ’.I 2 ./. 13 X 
it" NT= • ,F 4 . 0 f /. 13 X,*UR 0 ER= » . I 2 » / , 1 3 X , »COP= # ,I 2 .//) 

I F ( INPUT. EO. 2 ) GO TO 12 
IF(Y.LE.O.C) GO TO 6 
IF(YT .LE.Y.AN 0 .K.LT. 2 ) GO TO 8 
9 DENOM^Y T-Y T 2 

SL 0 PE 1 =( A 9 PK-ABPK 2 I/DEN 0 M 
SLOPE 2 = ( ABSK- AB SK 2 ) /DENOM 
SLQPE 3 =( ABNK- ABNK 2 ) /DENCM 
SLOPEA=( APHK-APHK 2 ) /DENOM 
SL 0 PE 5 =( APNK-APNK 2 ) /DENOM 
01 =ABPK-SLOPE 1 *YT 
82 =ABSK-SL 0 PE 2 *YT 
B 3 =A 6 NK- SLOPE 3 *YT 
84 =APHK-SL 0 PE 4 *YT 
b 3 =apnk -slope 5 *yt 
ABPT = SL 0 PE 1 *Y+B 1 
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TABLE IV-2. ( Coat'd) 


ABST=SL0PE2*Y+62 
A6NT = SLOPE2*Y+83 
APHT=SL0PE4*Y+B4 
APNT= SLOPE 5*Y +B5 
I F l INPUT. EO. 3 > GO TO 3 
GO TG 52 

6 kt ADI 5» 507 ) YT,ABPK,ABSK,ABNK ,APHK , APNK , VC I T 
C I**************** ******** A* ****** $*#**$***❖***>*><<★******** ****#**#$ 


C * READ IN TABULAR VALUES FOR Y-0.0 (NOT REQUIRED IF INPUT=2) * 
C * * 
C * A6PK IS THE BURNING AREA IN THE PORT IN IN**2 * 
C * ABSK IS THE BURNING AREA IN THE SLOTS IN IN**2 * 
C * A6NK IS THE BURNING AREA IN THE NOZZLE END IN I N**2 * 
C * A P HK IS THE PORT AREA AT THE HEAD END IN IN**2 * 
C * APNK IS THE PORT AREA AT THE NOZZLE END IN IN**2 * 
C * VC IT IS THE INITIAL VOLUME Of CHAMBER GASES ASSOCIATED WITH * 
C * TABULAR INPUT IN IN**3 * 


C **************** ** * ********* ********* ********************* ********** * 
507 FORMAT! 6X , F6. 2, 1 OX ,E 1 1 . 4 , 1 OX, El 1 . 4 , 8 X , E 1 1 .4 , / , 2 2X , E 1 1 .4 , <9X , El l . 4, 

1 8X » E 1 1 .4 ) 

WR I T E I 6 » 6 1 C > 

610 FORMAT! 13X , 'TABULAR VALUES FOR YT EQUAL ZERO READ IN') 

WRITE (6, 582) ASPK, ABSK ,ABNK ,APHK ,APNK. 

583 FORMAT! 13X, »ABPK*», IP E 11 . 4 , 5 X ♦ » AB SK= * , 1 PE 1 1 . 4 , 5 X , * ABNK= * ,i PEI 1.4, 

1 5X, , APHK = *, 1PE11.4,5X» * APNK= « , 1 PE i 1 .4 ) 

WRITE! 6 ,584 I VC I T 

584 FORMAT! 13X , *VC I T = * , 1PE11.4 ,// ) 

A8PT = A B P K 

ABST =AB SK 
ABNT = ABNK 
AP HT =AP HK 

apnt=apnk 

YT 2=YT 

IF! INPUT. EQ. 3 > GO TO 3 
GO TO 52 
8 YT2=YT 

A3PK2-ABPK 
A BNK 2 = A BNK 
AB SK 2=A B SK 
AP HK 2 - A PHK 
APNK 2 = A PNK 

HEAD! 5, 5051 Y T,ABPK ,AB SK ,ABNK ,APHK , APNK 
C ********************************************************************* 
C * READ IN TABULAR VALLES FOR Y=Y (NOT REQUIRED FOR INPUT=2) * 

C *** ****************************************************************** 
505 FORMAT! 6X , F6. 2, 10X , E 11 . 4, 1 OX , E 11 . 4 , 8X , £11 .4 , / , 22X, El 1 .4, 9X , Ell .4) 
WRITE (6,611) YT 

611 FORMAT! ///, 13X, 'TABULAR VALUES FOR YT= •»F7.3,» READ IN') 
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TABLE IV-2. (Cont'd) 


WKITE(6»583) A3PK , ABSK »ABNK » A PHK , APNK 
GU TO 9 
12 ABPT=0.0 
A i N T = 0 • 0 
ABST=0.0 

3 I F ( GR A I N ,NE ,2 ) GO TO 4 
AHPC=G.C 

ABNC = 0 . 0 
AHSOO.O 
GO TO 7 

4 IF(Y.LE.O.O) Rfc AO ( 5? 5 C l ) DO ,DI , DE LO I ,S » THETAG , LGCI t LGNI ♦ THETCN* THE 
ITCH 

£ fc*****^'-***** ***> (C# #*******#****$***♦*#* ************** ************* 

C * READ IN BASIC GEOMETRY FOR C.P. GRAIN (NOT REQUIRED FOR * 

C * STRAIGHT STAR GRAIN) * 

C * DO IS THE AVERAGE OUTSIDE INITIAL GRAIN DIAMETER IN INCHES * 

C * D I IS THE AVERAGE INITIAL INTERNAL GRAIN DIAMETER IN INCHES * 

C * DELDI IS THE DIFFERENCE BETWEEN THE INITIAL INTERNAL GRAIN * 

C * DIAMETER AT THE NOZZLE END AND 01 IN INCHES * 

C * S IS THE NUMBER OF FLAT BURNING SLOT SIDES (NOT INCLUDING * 

C * THE NOZZLE END) * 

C * THETAG IS THE ANGLE THE NOZZLE END OF THE GRAIN MAKES WITH * 

C * THE MOTOR AXIS IN RADIANS * 

C * LGCI IS THE INITIAL TOTAL LENGTH OF THE CIRCULAR PERFORATION * 

C * IN INCHES * 

C * LGNI IS THE INITIAL SLANT LENGTH OF THE BURNING CONICAL * 

C * GRAIN AT THE NOZZLE END IN INCHES * 

C * THETCN IS THE CONTRACTION ANGLE OF THE BONDED GRAIN IN RAD, * 

C * THETCH IS THE CONTRACTION ANGLE AT THE HEAD ENO IN RADIANS * 

Q **************** ***************************************************** 

501 FORMAT! 5X,F7. 3 , 6X , F7. 3 , 9X , F 7. 3 , 5X ,F 4. 0 » 9X , F7 . 5 , / ,7 X , F7 . 2 , 7X , F6 . 2 , 9 
lXt F7.5 t 9X,F7. 5) 

IF(Y.LE.O.O) WRITE (6, 601) 00,DI,DEL0I , S t THETAG , LGCI ,LGNI ,THETCN,TH 
1ETCH 

601 FORMAT ( 20X * *C ,P« GRAIN GEOMETRY* »/»l3X» *00= * ,F7.3 ,/ ,13X t * 01= ',F7 
1.3./.13X, •DELDI* • ,F7.3,/*13X,*S= • ,F4. 0 , / , 1 3X , • THET AG= •♦F7.5,/»l 
23X,*LGCI= *,F 7.2,/,13X,»LGNI= • » F6. 2 , / , L3X , * THETCN* • . F7 .5, /» 13X. • 
3THETCH= ■ v F7.5t//) 

D0SQ0=00*D0 

DISQD=DI*DI 

3NUM=ANUM*D0SQD 

TLL=TL 

I F ( OP-DtK • GE,3 ) TLL = 0.0 
YD I = 2 ,*Y+D I 
V D 1 $00= YD I* YD I 
ABSC=S*ANUM*{ OOSQD-YDISQD) 

IF( ABSC.LE.O.O) A8SC=0.0 
IF(YDI.GT.DO) GO TO 100 
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TABLE IT-2. (Cont'd) 


!F(THETAG.GT. 0,08727) GO TO 101 
IF(COP.EQ.C) GO TO 700 
IF(COP.EQ.l) GO TO 701 
I F ( COP • EQ . 2 ) GO TO 702 

LGOL GC I-( SQR T ( DOSQO-O I SOD ) - SORT { DOSQD- YD 1 SQD ) ) /2. -Y*C0T AN ( THcTCN ) 
GO TO 710 

702 LGOLGC l-( SQR T< OOSGD-D ISQD >-SGRT ( DOSQD-YD I SQD ) ) 

GO TO 710 

701 LGOLGC I-( SQRT{DOSQD-<OI+OELDI )**2)-SQRT(D0SQD-<YDI*DELDI )**2) )/ 

1 2.-Y*COTAN( THETCH) 

GO TO 710 

700 LGOLGC i-Y* (CO TAN ( THETCN J+COTAN ( THE TCH >) 

710 AI3PC«Pl*YDI*< LGC-TLL-S*Y) 

ABNC-0 . 0 
GO TO 732 
101 CONTINUE 

IF IC0P.EQ.0.0R.C0P.EQ.1) GO TO 720 

ABPC=PI*YDI*( LGCI-< SORT (DO SQD-D I SQD) -SORT (DOSQD-YD I SOD) )/ 

1 2.-TLL-< S+TAN( THE TAG/2 .)) *Y) 

GO TO 720 

7 20 ABPC=PI*YOI*< LGCI-Y*COTAN( THE TCH) -TLL- ( S+TAN(THETAG/2. )) *Y ) 

730 1 F ( COP * EQ . 1 .OR .COP . EG. 2) GO TO 731 
ABNC=Pl*(LGNI-Y*COTAN( THE TAG+THETCN) -Y*TAN (THETAG/2. ) )*( 01+ 

1 OELOI+Y+LGN I*SIN(THETAG) + Y*SIN(THETCN)/SIN<THETAG+THETCN) ) 

GO TO 722 

731 P.3 = l (DI+DELDI ) /2 . +LGN I * S I N ( THE TAG ) ) *C OS ( THE! AG ) -S I N ( TH£T AG ) * 

1 SORT ( <00/2. ) ^^2- ( (DI+DELOI ) /2.+LGNI*SIN( THETAG) )**2> 

A iNC = PI«(LGNI -R3-Y*TAN< THE TAG/2. ))*UDI + D6LDI) /2. + SQRT < < 00/ 

1 2. )< [ #2-(R3+Y)**2)*SIN( THETAG > +Y+ (R3+Y) *COS ( THETAG) ) 

732 IF(ARPC.LE.O.O) ABPC^O.O 
I F ( ABNC .L E .C. 0 ) ABNC=0.0 
GO TO 5 

100 ABNC- 0.0 
AbPC-0.0 

5 AP HT=ANUMY( DI +2.*RHT>**2 
IFIAPHT .GE.BNUM) APHT=BNUM 
IF<K.LT.2) APHT1=APHT 
APNT=ANUM*(CI +DEL0 I+2.*RNT ) **2 
I F < APNT. GE.BNUM) APNT = BNUM 
IF (GRAI N.NE.l ) GO TO 7 
ABP$=0.O . 

ABSS=0.0 
ABNS=G. 0 
GO TO 50 

7 IF(Y.LE.O.O) RE AO (5*5021 M S » LGSI » NP » RC ,F I LL » NN 

c ********** ****************************************************** ***** 

C * READ IN BASIC GEOMETRY FOR STAR GRAIN (NOT REQUIRED FOR * 

C * STRAIGHT C.P. GRAIN) « 


• 39 - 



age 17 


TABLE IV-2. (Cont'd) 

C * NS IS THE NUMBER OF FLAT BURNING SLOT SIDES (NOT INCLUDING * 


C * THE NOZZLE END) * 
C * LGSI IS THE INITIAL TOTAL LENGTH OF THE STAR SHAPED * 
C * PERFORATED GRAIN IN INCHES * 
C * NP IS THE NUMBER OF STAR POINTS * 
C * RC IS THE AVERAGE STAR GRAIN OUTSIDE RADIUS IN INCHES * 
C * FILL IS THE FILLET RADIUS IN INCHES * 
C * NN IS THE NUMBER OF STAR NOZZLE END BURNING SURFACES * 


C *#*******$ jjt*###*##****#***** *#*#$##*)(£$*:$(#$**£$***£****##!{<***>>**#*#*## 

502 FORMAT! 5 X , F 4. 0 , 7X , F7. 2 , 5X , F4. 0 , 5 X , F 7. 3 , 9X , F7 .3 , 3X , F4. 0 > 
IF(Y.LE.0«0) WR I TE { 6 ,602 > NS, LG SI »NP,RC»FILL,NN 
602 -FORMAT! 20X, 'BASIC STAR GEOMETRY* ,/, 13 X NS* • , F4. 0 , / ♦ 13X , • LGS I = *, 
1F7 . 2, / , 1 3X , *N P = *,F4.0,/, 13X,*RC= * »F7*3»/»13X,'FILL= *,F7.3,/,13X 
2, • NN= • , FA ,G, //) 

P I CNP =P I /NP 
RC SQD=RC*RC 
FY -F ILL +Y 
FY SQ0 = FY*FY 

IF ( STAR .FQ. I) GO TO 20 
I F ( STAR .EO. 2) GO TO 201 
IF(Y.GT.O.O) GO TO 179 

R E AD { 5, 42 1 ) TAUWW, L 1 »L2, ALPHA1 , ALPHA2 t HW 
C ***** ****** ***************** ************* *********** ****** *** ******** 


C * READ IN GEOMETRY FOR WAGON WHEEL (NOT REQUIRED FOR STANDARD * 
C * OR TRUNCATEO STAR GRAINS) * 
C * TAUWW IS THE THICKNESS OF THE PROPELLANT WEB IN INCHES * 
C * LI AND L2 ARE THE LENGTHS OF THE TWO PARALLEL SIDES OF THE * 
C * TWO SETS OF STAR POINTS IN INCHES * 
C * ALPHA! AND ALPHA2 ARE THE ANGLES BETWEEN THE SLANT SIDES OF * 
C * THE STAR POINTS CORRESPONDING TO LI AND L2 , RESPECTIVELY,* 
C * AND THE CENTER LINES OF THE POINTS IN RADIANS * 
C * HW IS HALF THE WIDTH OF THE STAR POINTS IN INCHES * 


C ***** **************** **** *************** **** ****** ******** **** ******* 

421 FORMAT< 3( 6X,F5.2> ,2(10X,F7.5) ,6X,F5.2) 

WRITE! 6, 422) TAUWW, L i , L2 » ALPHA1 , ALPHA2 ,HW 

422 FORMAT( 20X , 'WAGON WHEEL GEOMETRY* ,/ ,13X ,« TAUWW* ' ,F5«2,/ , L3X, 

1 ’Ll* *,F3.2,/,13X,*L2= • ,F5.2 , / ,13X, • ALPHA1* » ,F7 .5,/ , 13X , 

2 * AL PHA2* *,F7.5,/,13X,*HW= ',F5.2,//> 

AL P2 = AL P FA 2 

XL 2*1 2 

L F W=RC- T AU'w W- F I LL 
LFWS«0=LFW*LF W 
THLTFW=ARSIN( ( HW+F I LL) /LF W 1 
SLFW=LFW*SIN( THETFW) 

179 KKK=0 
S G=0 • 0 

ENUM={RCSQD“LFWSQD-FYSQO) /(2.*LFW*FY) 

AL PHA2* ALP2 


-4d- 


VGE 18 


TABLE IV-2. (Coat'd) 


L 2=XI_ 2 

190 Y TAN = Y* TAN < ALPHA2/2.) 

CO SAL P = COS { ALPHA 2 ) 

S INALP=SIN( ALPHA2) 

IF(YTAN.GT.L2 ) GO TO 182 
1 F < FY .G T . SL fW ) GO TO 181 

SGW=NP*(L2-2.*YTAN-MSLFW-FILL) /Si NALP-Y*COTAN( ALPHA2)+FY* 

1 (PID2+THETF W>+<LFW+FY)#( PIDNP-THETFW) ) 

GO TO 183 

181 IF(Y.GT.TAUWW) GO TO 184 
SGW=NP*{FY*(PIDNP+ARSIN( SLF W/F Y) ) + ( P I DNP-THETF W ) *L F W ) 

GO TO 183 

164 SGW«NP*FY*<THETFW+ARSIN( SLFW/FY)-ARCOS(ENUM) ) 

GO TO 183 

182 YPO=-SLFW 

IF( ALPHA2.G6.P (02) GO TO 222 
Q--FILL+L2*TAN< ALPHA21-Y/COSALP 

X P l = ( -Q*T AN( AL PHA2 ) -SQRTT -Q*Q+FYSQD/COSALP*CQ$ALP) ) *C0SALP*C0SALP 
YP I=XPI*TAN<ALPHA2)+Q 
XPO= ( YPO-Q ) *CQT AN ( ALPHA2) 

GO TO 223 

222 XP I =Y-L 2 

YP I a- SORT ( FYSOD-XPI*XPI) 

XPG = XP I 

223 FYL$=SQRT( SLF W* SL F W+XP I * XPI ) 

XP I02=( XP I -X 3 0 ) * ( XPI-XPO) 

YP I0 2={ YP I-YP0)*{ YPI-YPO) 

I FI FY.GT.FYLS ) GO TO 186 
IF(Y.GE.TALUW) GO TO 185 

SGirt»NP*( SORT( XPI02+YPI02)+FY*(PID2+THeTFW-AR$IN (XP I/FY) )+( LFW+FY )* 
1 (P IDNP-THETFW) > 

GO TO 183 

185 SGw=NP*( SQRTl XP 1 02 *YP 1 02 ) +FY* ( P I D2-ARSI N ( XP I /FY ) -ARCOS ( ENUM >) ) 

GO TG 183 

186 IF(Y.GT .TALK*) GO TO 187 
SGW=NP*(FY*(PIDNP+ARSIN< SLFW/FYJ ) + (PI DNP-THETF W)*LFW) 

GO TO 183 

187 SGU=NP*FY*(THETFW+ARSIN( SLFW/FY)-ARCOS(ENUM) ) 

183 IF(SGW.LE.C.O) SGW=0.0 
IF(Y.GT.O.C) GO TO 188 

AGS2-«5*(PI*RC SQD-NP* LFW#SLFW ! MCQS<TH£TFW)-*SIN(THETFW)*CQT AN (ALPHA 

1 2)-2.*(L2+FlLL*TAN(ALPHA2/2. ) ) /LFW) -< Pl-THETFW*NP) *LFWSQ0-2.*NP*F 

2 ILL*(L2 + SLFW/SINALP + LFW*(PIDNP-THETFW) + (PIDNP+PID2-1./SINALP> * 

1 FILL/2.)) 

AG$= AGS+AGS2 

188 CONTINUE 

S G = SG + SGW 

IF(KKK.EO.l) GO TO 24 


-4l- 



TABLE IV-2. (Cont'd) 


L2=L 1 

AL PHA2= ALPhAl 
KKK = 1 
GO TO 190 

201 IF(Y.LE.O.O) RFAD(5t503) RP,TAUS 

£ ***** ***** ** * *** ***************************************************** 


C * READ IN GEOMETRY FOR TRUNCATED STAR (NOT REQUIRED FOR * 
C * STANDARD STAR OR WAGON WHEEL) * 
C * RP IS THE INITIAL RADIUS OF THE TRUNCATION IN INCHES * 
C * TAUS IS THE THICKNESS OF THE PROPELLANT WEB AT THE BOTTOM * 
C * OF THE SLOTS IN INCHES * 


C * **** *** ****** ** ***** ******* ******** ********** **** **************** *** 
503 FORMAT! 5X, F7. 3, 7X,F7.3) 

IF(Y.LE.O.O) WR I TE ( 6» 603 ) RP,TAUS 

603 FORM AT ( 20 X » * TR UNCATE 0 STAR GEOMETRY* ,/ fl3X t * RP= « » F7. 3, / , 13X f * T AUS 
1= *,F7.3t//> 

THETA$=PIONP 

RPY-RP+Y 

LS=RC-TAUS-FILL-RP 
KPL=RP + LS 

THETSl=THETAS-ARSIN(FY/RPY) 

IF(THETSl.LE.O.O) GO TO 110 
IFtY.LF .TAOS) GO TO 103 

T HFT AC- AR S IN( ( RC SQD-R PL*RP L-F YSQD ) / < 2 . *F Y*RPU) 

IF(THETAC.Gfc.O.O) GO TO 104 
IF< Y.LT.RC-RP ) GO TO 105 
S G=0 . 0 
GfJ TO 1 A 

103 SG=2 «*N P*( R PY*THET Sl+L S- ( RPY*COS (THE TAS-THETSl ) -RP)+PID2*FY) 

GO TO 14 

104 SG = 2.*NP*(RPY*THETS1+LS-(RPY*C0S< THETAS-THETS1 ) ~RPH-FY*THETAC) 

GO TO 14 

105 SG-2. *N P*(RPY* THE TS 1+SQRT ( RC SQD-F YSQD) -SORT ( RPY*RPY-FYS QO ) ) 

14 IF(Y.LE .0.0) AGS=PlMRCSQO-RP*RP)-NP*(PI*FILL*FILL/2.+2.*LS*FlLL) 

GO TO 31 

110 THETAF=THETAS 

THETAP=2.*THETAS 
TAUWS=T ALS 
GO TO 111 

20 IF(Y.LE.O.O) READ ( 5» 504) THE TAF , THE TAP,TAUWS 

C ***X< ****** ****** ***************************************************** 


C * RE AO IN GEOMETRY FOR STANDARD STAR (NOT REQUIRED FOR * 
C * TRUNCATED STAR OR WAGON WHEEL) * 
C * THETAF IS THE ANGLE LOCATION OF THE FILLET CENTER IN RADIANS * 
C * THETAP IS THE ANGLE UF THE STAR POINT IN RADIANS * 
C * TAUWS IS THE WEB THICKNESS OF THE GRAIN IN INCHES * 


C ********************************************************************* 
504 FORMAT! 9X, F7. 5» 9X,F7.5,8X t F6.3> 



TABLE IV-2. (Cont'd) 


IF(Y.LE.O.O) WR I TE < 6» 604 ) THE TAF , THE TAP , TAUWS 
604 FORMA T( 20X» 'STANDARD STAR GEOMETRY* ,/» 13X» ' THETAF= ' t F7 • 5* / » 13X . 1 T 
1HE TAP= • , F7.5 ,/, 13X, *TAUWS= »,F6.3»//) 

THETAS=PI/NP 
THETS 1= 1 .00 
111 LF=KC-TAUWS-F ILL 
CNUM = < Y+F ILL ) /LF 
DNUM=SIN( TFETAF > /SIN (THE TAP/2. ) 

FNUM= ( RC SOD-L F*LF-FYSQD)/(2.*tF*FY) 

FNUM=SIN( THETAF) /COS( THE TAP/2.) 
i F ( CHUM . LE .FM UM ) GO TO 106 
I F ( Y . Lfc . I AUViS ) GO TO 107 

SG=2.*NP*FY*( THETAF+ARSINt SIN(THETAF) /C NUM) -ARC OS < ENUM ) ) 

GO TO 23 

106 IF(Y.LE.T A UWS ) SG=2 . *NP* L F* I DNUM+CNUM* ( PI 02+THE T AS-T HE T AP/ 2 . 
1-COTAN ( THFTAP/2.) M-THETAS-THETAF ) 

IF(Y.LE.TAUUS) GO TO 23 

SG=2 . *NP*{ FY<* AR S IN ( ENUM- (THETA S-THETAP/2.) ) +LF*DNUM-FY*COTAN( THETA 
IP / 2. ) ) 

GO TO 23 

107 SG=2.*NP*LF*< CNUM* ( THE TA S+AR S I N < SI N (THE TAF ) /CNUM) ) +THETAS-THETAF) 

23 IF (THETSL.LE.O.O) GO TO 14 

IF(Y.LE.O.O) AGS«PI*RC*Q2-NP#LF*LF*<SIN(THET.AF)*(C0S(THETAF)- 
1S IN( THE T AF )*C QIAN ( THE TAP /2 . ) ) + THE TA S-THE TAF+2 . *FI LL/LF* ( S I N ( THET AF 
2)/SlN(THFTAP/2.)+THETAS-THETAF+FILL/(2.*LF)*(PI D2+THET AS -THE 
3TAP/2.-CUTANI THETAP/2.) ) ) ) 

24 CONTINUE 

31 IF(SG.LE.O.O) SG=0.0 

IF(K.EQ.O.OR.K.EQ.2) SGN=SG 

IF(K.LE.l) sgh-sg 

IF(Y.LE.O.C) SG2=SG 

I F { K . FO . 2 ) GO TO 37 

RAVEDT=R1+(SG+SG2)/2.*RBAR*DELTAT 

RNDT=R2+( SG+SG2)/2.*RNAVE*DElTAT 

RHDT=R3+($G+SG2)/2.*RHAVE*DELTAT 

R1=RAVEDT 

R2=RNDT 

R3=KFDT 

SG2-SG 

GO TO 3 8 

37 IFtKOUNT.NE.l ) GO TO 39 
SG3=SG 
R 4=R 1 
R5=R2 
R6 *R'J 

39 RAVEDT=R 4+ ( SG +SG3 ) / 2. *RB AR *DE LTA T 
RNDT=R5+{ SG+SG3)/2.*RNAVE*DELTAT 
RHDT=R6+( SG+SG3 )/2.*RHAVE*DELTAT 
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RHDT -R6 + ( SG+SG3)/2.*RHAVE*DElTAT 

P.4=R AVFOT 
R 5=RN0T 
R 6=RHDT 
SG3=SG 

38 ABSS=(AGS-RAVFOT)#NS 

IF(ABSS.LF .O.O.OR.SG.LE.O.O) ABSS=0.0 
ABNS=( AGS-RNDT )*NN 

IF! AbNS.LF .0. O.OR.SG.LE.O.O) ABNS=0.0 
IFiUROER.LE.2) ABPS=(IGSI~Y*(NS+NN) 1*SG 
I F { 0RDLR .LE .2 ) GO TO 36 
A8PS=(LG$I-Tt - YMNS+NN ) ) * SG 
36 PIRCRC=P I*RCSQD 

AP N$-P l RCRC-A GS+RHDT 

IFCAPHS.GE.PIRCRC.OR. SG.LE.O.O) aphs^pircrc 
APNS = P I RCRC-A GS+RNDT 
IFtK.LT. 2) APHS1=APHS 
IF(APNS.GF.PIRCRC) APNS=PIRCRC 
50 IF(NT.EQ.O.O) GO TO 37.1 

IF!Y.LE.C.G) REAO ( 5» 5C 6) LTP,DTP ,THE TTP tTAUEFF 
C, **************** ************ **** ***************************** ******** 


C * REAO IN GEOMETRY ASSOCIATED WITH TERMINATION PORTS (NOT * 
C * REQUIRED IF NT*OJ * 
C * L TP IS THE INITIAL LENGTH OF THE TERMINATION PASSAGES * 
C * IN INCHES * 
C * DTP IS THE INITIAL DIAMETER OF THE TERMINATION PASSAGE * 
C * IN INCHES * 
C * THETTP IS THE ACUTE ANGLE BETWEEN THE AXIS OF THE PASSAGE # 
C * ANC THE MOTOR AXIS IN RADIANS * 
C * TAUEFF IS THE ESTIMATED EFFECTIVE WEB THICKNESS AT THE * 
C * TERMINATION PORT IN INCHES * 


C ***** *********** ***************************************************** 
506 FORMAT! 7X,F6. 2 » 7X , F5. 2 ♦ 1 OX ,F 7. 5 , 1 OX *F6. 3 1 

IF(Y.LE.O.O) WR I TE ( 6» 606 ) LT P , DTP , THETTP , TAUEFF 
606 FORMAT! 20X, 'TERMINATION PORT GEOMETRY' ,/ t 13X ,* LTP= * , F6- 2t / ♦ 13X ♦ ' D 
1 TP * ' »F5.2»/» 13X» 'THETTP® • , F7. 5 , / , I3X, • TAUEFF= » t F6.3 ? //> 

DAtST =NT*3. 141 59* { ( OTP+2. *Y> * I LTP-Y/SI N ( THETTP l ) -(DTP+2. *Y) **2/4 .+ 
HY+OTP/2. >* (DTP/2. ) *11 .-1«/SIN( THETTP) ) ) 

IF(Y.GE. TAUEFF) DABT=0.0 
371 IF(Y.GT.O.O) GO TO 52 
I F < NT .NE.O.C) GO TO 45 
LTP=0.0 
DT P= 0 . 0 

45 IF! GRAIN .NE.2 ) GO TO 49 
LGC1 = 0» 0 
LGNl-0.0 

49 I F ( GRAIN. EG, 1 ) IGSI*0.0 

VC 1*1. 1*( ANUM*DISQD*(LGCI+IGNI ) + ( ANUM*DOSQD-AGS ) *LGS I *NT *LT P*ANUM* 
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1 QTP*DTP)+VC IT 

52 BBP=0.0 
03 S=0 . 0 
BRN=0.0 

ABPGRT=A8PT+A8PC +ABP S+DAB T+BBP 
AOSLOT = ABST+ABSC + ABSS'*'BBS 
ABNOZ=ABNT ♦ AB NC +ABNS+B BN 
SUHAB=AOPGRT+ABSLO T+ABNOZ 
IF(K.EQ.O) GO TO 99 

IRK .EG. 1 ) ABMAIN=ABPORT+ABSLOT+ABNOZ 
K = K + L 

IF(K.GT.2> GO TO 69 

GO TO 2 

69 ABTO=ABPORT+A BSLOT+ABNOZ 
99 CONTINUE 

IF(Y.GT.C.C) GO TO 70 
ABP i = ARPORT 
ABN I = ABNOZ 
AOS 1 = AB SLOT 

70 AbP2-(ABPi+ABPQRT)/2. 

ABN2=( ABNI+ABNOZ) /2. 

AB$2 = (ABSl+ABSLQT)/2* 

I F { INPUT. EQ. I > GO TO 76 
GO TO ( 71,72, 73, 74), ORDER 

71 APHEAD^APHSl 
APNOZ=APNT 
SG=SGH 

GO TO 75 

72 AP HEAD= A PHT 1 
APNOZ =A PNT 
SG=0 . 0 

I F ( GR A I N . EO .3 ) SG= ( SGH + SGN ) 11 . 

GO TO 75 

73 AP HE AD = APHT 1 
APNOZ =A PN S 
SG*SGN 

GO TO 75 

74 APHE A0-APHS1 
APNOZ -APNS 
SG=SGN 

GO TO 75 
76 AP HEAD- APHT 
APNOZ -APN T 

75 Y=YB 

01 FF=SUMAB~SUM2 
DADY = DI FF/CEL Y 
ABP 1=AGP0RT 
ABN1=ABN0Z 
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AR S l=AB SLOT 

IF(ZW.GE.O.O) GO TO 77 
A3MUABMAIN 
ABMA IN=A6TG 
ABTG=ABMl 
77 RETURN 
END 


(Cont r d) 
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SUBROUTINE OUTPUT 

C ** ****** **< ;*** ** ***** ***************** ************ ******************* 


t * SUBROUTINE OUTPUT CALCULATES BASIC PERFORMANCE PARAMETERS * 
C * AND PRINTS THEM OUT * 
C * (WEIGHT CALCULATIONS ARE PERFORMED IN THE MAIN PROGRAM) * 
C * T IS THE TIME IN SECS * 
C- * V IS THE DISTANCE BURNED IN INCHES * 
C * RNOZ IS THE NOZZLE END BURNING RATE IN INCHES/SEC * 
C * RHEAD IS THE HEAD END BURNING RATE IN INCHtS/SEC * 
C * PONOZ IS THE STAGNATION PRESSURE AT THE NOZZLE END IN PSIA * 
C * PHEAD IS THE PRESSURE AT THE HEAO END OF THE GRAIN IN PSIA * 
C * PTAR IS THE PORT TO THROAT AREA RATIO * 
C * MNOZ IS THE MACH NUMBER AT THE NOZZLE END OF THE GRAIN * 
C * PATM IS THE ATMOSPHERIC PRESSURE AT ALTITUDE IN PSIA * 
C * SUMAB IS THE TOTAL BURNING AREA OF PROPELLANT IN IN**2 * 
C * SG IS THE BURNING PERIMETER IN INCHES OF THE STAR SEGMENT * 
C * ( IF ANY) * 
C * CFVAC IS THE VACUUM THRUST COEFFICIENT * 
C * FV AC IS THE VACUUM THRUST IN LBS * 
C * F IS THE THRUST IN LBS * 


C ********** ******************************** *>)(** *************** ******** 
REAL MGEN,MDI S,MNOZ,MNl, JROCK , N , L , MEl , ME , I S P , IT OT , MU , MASS , ISPVAC 
REAL M2,MCBAR, ISP2, ITVAC 
COMMON /CONST 1/ZW,AE, AT, THE TA,ALF AN 
C0MM0N/C0NST2/CAPGAM,ME,BQT,ZETAFtT8»HB»GAM 

COMMON/ VAR I AI/Y,T,DELY, DELTA T, PONOZ, PHEAD, RNOZ, RHEAD, SUMA8,PHM AX 
COMMON/ VAR I A3 / I TOT, IT VAC , JROCK , I SP , I SPVAC ,MDI S , MNOZ , SG , SUMMT 
COMMON/ VAR IA5/ABMA IN, ABTO,SUMDY, VC I , VC 
CUMMGN/PLOTT/NUMPLT( 16) , I PO*NDUM, NP » I OP 

DIMENSION TPL 0 T ( 200 ) »PNPL0T(2 00) ,PHPLOT(200) ,FPL0T(200) , 

1 FVPLOT ( 200) tPNPLQT (200) ,RHPL0T(20G> ,YBPL0T(200) , ABPLOT ( 200) , 

1 SGPLOT< 200) , VCPLOTI200) 

DATA G/ 32 • 172 5/ 

I F ( NUUM.EO • 1 ) GO TO 2 

NP^NJP + I 

Y 3 =Y 

iF(Y.LF.O.O) M2=MDIS 
MDBAR = ( M2+M0I S)/2. 

SUMMT =SUMM T+MOBAR*DELTAT 
PTAR = 1. /JROCK 

PRES= ( i.+eCT/2.*ME*HE)**(-GAM/R0T) 

ALT=HB*( T / TH) *# ( 7./3. ) 

PATM=L4.69£/EXP (0.431 03E-04* ALT) 

IF C MO IS • LE. 0.0. OR. PONOZ. LE. 0.0) GO TO 45 

CF=CAPGAM*SCRT( 2.*GAM/BOT*(1.-PRES**(BOT/GAM) ) )+AE/AT*( PRES-PATM/P 
10M0Z ) 

CFVAC=CF+AE/AT*PATM/PONOZ 

F=Z ET AF*CO S( THE T A )* PONOZ* AT* ( ( 1 . +COS { ALF AN) ) /2. *CF+ ( l.-COS ( ALFAN) > 
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1/2»*AE/AT*(PRE S-PATM/PONQZ) ) 

IF(F.LE.O.O) F = 0.0 
IF { Y «LE ,0.0 ) F2=F 
FBAR = < F + F2I/2. 

FVAC=ZETAF*CO S( THETA)*P0N0Z*AT*((l.+C0S(ALFAN))/2. *CFVAC+ ( 1 .-COS( A 
1L FAN ))/2.*AF/AT*PRES) 

IF(Y.LF.O.O) FV2=F VAC 
FV BAR = ( FV2+FVAO/2. 
lSP=F/( MDIS*G> 

I S P V AC= FVAC /I MDIS*G) 

lTOT=ITOT+FBAR*DELTAT 

ITVAC=ITVAC*FVBAR*DELTAT 

M2=MDIS 

F2 = F 

FV 2=F VA C 

IF(PHEAO.GT.PHMAX) PHMAX=PHEAD 
GO TO A 7 
45 CFVAC=0.0 
FV AC = 0. 0 
F = G.O 

47 WR ITE<6, L) T, Y0 , RNOZ , RHE AO , PONOZ , PHE AO » PTAR , MNOZ,S UMA8 , S G, PATM t CFV 
14C » FVAC » F 

1 FORMAT* //, 13X ,♦ TIME* « ,F7 . 3 , 1 2 X, • Y= » ,F6.3 ,/ ,13 X,' RNOZ= »,1PE11.4, 
1* RHEAD= »,1PE11.4,» PQNQZ= ',1PEU.4,' PHEAD= • , 1 PEI 1.4, / , 1 3X , 
2' PTAR= ' , 1PEI 1.4, * MNOZ* »,lPE11.4t* SUMAB= ',lPE11.4, f SG= 

3 • , LPE1 1.4,/, L3X,'PATM= ',1PE11.4,' CFVAC= ',1PE11.4,» FVAC= ', 
41PE11.4,' F= SLPE11.4) 

► IF(IPQ.FQ.C) RETURN 

► TPL'DT ( N P ) = T 

► PNPLOTI NP )= PONOZ 

► PHPLOTI NP )=PHtAD 

► F P L ‘ ) T ( N P ) = F 

► FVPLGTt NP)=FVAC 

► RNPLOT(NP)=RNOZ 

► RHPLOTC NP )=RHEAO 

► Y B PLOT ( NP i =Y8 

► ABPLOT( NP)=SUMAB 

► SGPLUT* NP > = SG 

► VC PLOT { NP )=VC 

► RETURN 

► 2 NP=NP+2 

► I0P = 1 

► no 10 04 1=1,16 

► I F l NUMPL T ( I ) • EQ • 1 ) GO TQ 1003 

► GU TO 1004 

► 1003 GO TO ( 10, 20, 30,40,50,55,60,70,75,80,90,95,97,100, 110,115) , I 

► 10 CALL PLOTITITPLOT, 'TIME I SEC S ) ' , 1 1 , PHPLOT , • PHFAO <PSJA)',12» 

► 1 PNPLOT, 'PONOZ* »5,NP»l»*DUMMY*,5) 
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> GO TO 1004 

> 20 CALL PL GT I T( T PLG T» ' T I ME { SEC S ) • , LI , PNPLOT, ' PQNOZ ( PS I A ) * , 12, PHPLOT 

> i f \PHEAD ( PSIA l ' , 12,NP,1, 'DUMMY* ,5) 

> GO TO 100 A 

> 30 CALL PLOTI T( TPLOT, 'TIME ( SEC S ) ’ , 1 1 » PH PLOT , * PHE AO* , 5 , PNPL OT 

> 1, ' PONOZ ', 5, NP ,3, ’PRESSURE (PSIA)*, 15) 

> GO TO 1C04 

> 40 CALL PLOTIT(TPLOT, 'TIME (SECS) ' » 1 1 , RHP LOT , 'RHEAO (IN PER SEC)M8» 

> 1 PH PLOT » ' PHE AD ( P SI A ) ' , 12tNP *1 , •DUMMY 1 , 5) 

> GO TO 1004 

t> 50 CALL PLOTITt TPLOT, 'TIME ( SEC S ) • , 1 1 , RNPl OT , • RNOZ (IN PER SEC)*, 17, 

> 1PN PLOT , * PONOZ (PSIA) ' » 12 ,NP*» 1 » 'DUMMY* ,5) 

t> GO TO 1004 

> 55 CALL PLOTI T( TPLOT, * TIME ( SEC S) * » 1 1 » RHPLGT , ‘RHEAO* »5,RNPL0T , 

> 1 'RNOZ' ,4,NP, 3, 'BURNING RATE (IN PER SEC) ',25) 

> GO TO ICO 4 

> 60 CALL PLOTITtTPLOT, 'TIME ( SEC S ) ' , 1 1 , AB PLOT , ' TOT A L BURNING AREA (SQ 

> 1 IN ) ' , 26, PNPLOT , 'PONOZ' ,5»NP»l,' DUMMY* ,5) 

> GU TO 1CC4 

> 70 CALL PLOTITUPLOT, 'TIME (SECS) ', 11, SG PLOT, 'STAR PERIMETER (IN)', 19 

> 1, PNPLOT, *P0N0Z',5,NP,1,'0UMMY»,5) 

> GO TO 1CC4 

> 75 CALL PLOTI TITPLOT, ' TIME ( SEC S ) ' , 1 1 ♦ ABPLOT , ' TOT A L BURNING AREA (SQ 

> 1 IN ) ' , 26 » SGPLO T, 'STAR PERIMETER ( I N) • , 1 9 , NP ,2 , ' DUMMY ' ,5 ) 

> GO TO 1004 

> 80 CALL PLOT IT(TPLOT, 'TIME ( SECS) ' ,11 ,FPLOT ,' THRUST ( LBS )', 12 , PNPLOT , 

> 1 • PONOZ * ,5.NP» 1, 'DUMMY* ,5) 

> GO TO 1004 

> 90 CALL PLOTIT(TPLOT, ‘TIME ( SEC S 1 1 , F VPLOT VACUUM THRUST (LBS)', 19 

> 1, PNPLOT, 'PONOZ* , 5,NP, 1 , 'DUMMY' ,5) 

> GO TO 1004 

> 95 CALL PLOTIT(TPLOT,'TIME ( S EC S )', 1 1 , F PLOT ,' THRUST *, 6 , FVPLOT , 

> l 'VACUUM THRUST* , 13, NP, 3 , 'THRUST (L8S)',12) 

> GO TO 1 C04 

> 97 CALL PLOTITUPLOT, 'TIME ( SEC S )', 1 1 , VCPLOT ,' CHAMBER VOLUME (IN**3)' 

> 1 , 22, PNPLOT, • PONOZ' ,5,NP, 1, 'OUMMY* ,5) 

> GU TO 1004 

> 100 CALL PLOTIT(YBPLOT, '8URNED DISTANCE ( I N) • , 20 , A8 PLOT , ' TOT AL BURNING 

> I AREA (SO IN) ',26, PNPLOT, 'PONOZ' ,5, NP, l, 'DUMMY' ,5) 

> GO TO 1C04 

> 110 CALL PLOT I TIYBPLGT, 'BURNED DISTANCE ( I N ) ' , 20, SGPLOT , ' STAR PER l MET E 

> 1R UN )', 19, PNPLOT, 'PONOZ* ,5, NP,1 ♦•DUMMY* ,5) 

> GO TO 1004 

115 CALL PLOTIT(YBPLOT, 'BURNED DISTANCE ( I N ) ' , 20 , ABPLOT , • TOT AL BURNING 

> I AREA (SQ IN) ',26, SGPLOT, 'STAR PERIMETER ( I N) ' , 19, NP,2 , * DUMMY* , 5 ) 

> 1004 CGNT INUE 

RETURN 

END 
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SUBROUTINE IGNITN 

c * **** ******* **** ***************************************************** 


C * SUBROUTINE IGNITN CALCULATES THE PRESSURE RISE DURING * 
C * THE IGNITION PERIOD * 
C * AS IG IS THE IGNITER THROAT AREA IN IN**2 * 
C * WIGTQT IS THE TOTAL WEIGHT OF THE IGNITER IN LBS * 
C * MIGAV IS THE IGNITER AVERAGE MASS FLOW RATE IN SL/SEC * 
C * PC I G IS THE IGNITER PRESSURE IN IBS/IN**2 * 


COMMON/CON STi/ZW,AE , AT , THE TA , ALFAN 
COMMON /CUN ST2 /GAMMA, ME ,80T , ZE TAF , TB , HB , GAM 

COMMON/ VAR I Al / Y , T I G tOELY , DELTA T,PQNO*,PHE AD, RNOZ,RHE AD, SUMAB.PHMAX 
COMMON/ V AR I A2/ABP0RT,ABSL0T, ABNOZ, A PH EAD , APNGZ , OADY , ABP2 , ABN2, ABS2 
CCMMON/VAR1A3/ITOT, IT VAC , JROC K, I SP, I SPVAC , MOI S , MNOZ , SG,SUMMT 
COMMON/VARIA5/ABMAlN,ABTO,SUMDY, VC I ,VC 
COMMON/ I GN 1/<A,KB,UFS,RH0,A,XN,L, PMIG,TI 1,TI2,CS1G 
COMMON/ IGN 2 /ALPHA, beta ,pbig,rrig,deltig,x,top,zap 
► COMMON/ PLDTT/ NUMPL T ( 16 ) , IPO , NDUM , I P T, I OP 

DIMENSION 0(9 > 

PFAL K I 4 ) , L ,KA,KB, JROCK, J2,MIG,MI G4V , MSRM, Mfc , MD I S , MNOZ , MNOZ I,MNL 

DATA Ai»A2, A3, A4/,i7476»- .551481 ,1.205536 , .171185/ 

DATA ft { l>,e(2),B(3),B(4)*B(5)/0. , . 4 ,. 455737 , 1 , , .296978/ 

DATA ft( 6) , B< 7 ) ,B ( S ) ,B ( 9 > / . 1 58 76 , . 21 31 , -3 . 050965 , 3 . 8 328 64/ 

c ***** ********* ** ***************************************************** 
C * THE A 1 S AND Q*S ARE CONSTANTS FOR THE RUNGE-KUTTA INTEGRATION * 

C ************** ** ********************************** ** **************** * 
DATA G/32 • 1725/ 

X X . C5*P0N0Z 
I PLUG- 0 
PONOZ I = PONCZ 
RHEADI*RHEAD 
RNOZ I =RNOZ 
PHEAOl=PHEAD 
DfcLTT =DELT AT 
DISM=MDIS 
DELTAT-OELTir, 

SUMAB I = SUMAH 
MNOZI=MNOZ 
MNOZ -0.0 
PHEAD-0 .0 
« N 0 ?. = 0 . 0 
MO I S-0. 0 
AB 1 = 0 • 0 
T I G I = C . 0 
PC 1=14.696 

ug=o.o 

PCNEW=1 4 .656 
SUMAB=0 . 0 
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PC IG=14.696 
PHE A0=14 .696 
P0NGZ= 14.696 
SLOPE=SUMABI/L 
G2=GAMMA*GAMMA 
J 2=JR0CK* JROCK 
GJ=G2*J2/2. 

MIGAV=. 2*AT/G 

AS IG = 4. IGAV^C S IG/< 4.*PMIG— RRIG* (TI2-TI 1 ) ) 

W I GTOT = G*M I GA V*( 5. * < T I 2- TI l ) /6. ) 

WP, ITE(6i999 ) ASIG»WIGTOT»MIGAV 

999 FQRMAT( /// »2QX» * IGNITER SIZE CALCIJLATI ONS * » / » 13 X» * AS I G= * » F6.2* / ♦ 
l 13X» , WIGT0T= , »F7.2»/ »13X»*MI GA V= * ,F8.3t///J 
WR ITE(6, 10) 

10 FORMAT ( 3 3X , •**#***#****** I <‘**#* 1 « t *<‘*** J i' ; ** , ! t '»/ »33X * 1 **** IGNITION TRA 
1NSIENT **** « f / t 33X, =0= ^ ^ * * ) 

18 NNN*0 

CALL OUTPUT 
WR I T E l 6 » 30) PCIG 
30 FORMAT! 13X»'PCIG= *,1PE11.4) 

9 CONTINUE 
DO 8 N= 1 1 4 
IF(N.EQ.l) PC=PCI 
IF(N.E0.2) PC =PC I +8( 2 ) *K ( 1 ) 

I F ( N .EG . 3) PC =PC 1 + 0 < 5>*K, ( 1 ) + ft ( 6) *K(2) 

IFIN.E0.4) PC=PCI+8( 7)*K(1 )+8 < 8 J *K ( 2) +B ( 9 1 *K ( 3 > 

TIG=TIGI+S(N>*DELTIG 

SUMAB=ABI+SL3PE*UFS*B<N>*D£LTIG 

I F ( SUMAR.GT.SUMABI ) SUMA8=SUMABI 

PHEAD=PC 

IF (MO IS .NE .0. O') PHEAD=PC*( l.+GJ) 

REf. AC=A*PHEAO**XN 

IF(TIG.LE.TII) PCIG=PM IG*TIG/TI 1 

IFITIG.GT.TU .AND. PCIG. GT.PHEAD) PCIG=PMI G-RRI G*( T IG-T 1 1 ) 

IF( PC IG.LE .PHEAD) PCIG=PHEAD 
M I G-0 .0 

IF(PC IG. GT.PHEAD. AND. TIG. LE.TI 2/2. ) MIG-PCIG*ASIG/CSIG 

CS rR=KA+K8*PC 

MOIS=PC#AT/CSTR 

IF (PC .LE.PB IG .AND. IPLUG.EQ.O) GO TO 7 
I PLUG= 1 
MNQ7=MNOZ I 
PN OZ =PC*( l.-GJ ) 

L IT=MDI S*X/APNOZ 
RN 1=RHEAD 
AZ=ALPHA*Z IT**.8 
XL =U F S*T IG 
IF(XL.GT.L) XL=L 
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4 EX -XL** .2* EXP ( BETA*RN i*RH0/Z IT) 

KKO? ®RN l-< RNl-A*PNOZ**XN-AZ/EX) / (l.+AZ*BETA*RHO/(Z IT*£X) ) 
IFUBSIRNl-RNOZI.LE. 0.002) GO TO 5 
RNl®RNOZ 
GO TO 4 

7 M D IS = 0 • 0 
NNOZ=0.0 
PN ;)Z -PC 
RNCIZ=RHEAD 

5 CONTINUE 

KSRM=RHO*SUMAfi* (RNOZ+RHEAD) /2. 

DENOK = ( VC I / ( l 2.*CSTR*C STR*G2 ) ) * (1 . - (2 . *KB*PC ) /CSTR ) 

DPDT = (MIG+l*SRM-MDIS)/DENOM 

IF (f)PQT .LT.O. O.AND.PC.LT.20.0) 0P0T=0.0 

K ( N )=OELT IG*DPDT 

8 CONTINUE 

PCNEW=PC I+A1*K( l)+A2*K (2)+-A3*KI 3>+A4*K(4) 

PH EAU=PCNE W 

IF(MDIS.GT.O.O) PHEAD*PCNEW* < l.+G J ) 

PGNOZ=PCNEV. 

XX Y = AGS ( PONGZ- PQNOZ I ) 

IF ( PCNEW »L E .1 *001* PCI .AND. SUMAB. EQ. SUMABI .AND.XXY.LE.XXX) GO TO 13 

AB 1 = SUV, A 8 
T I G I® T I G 
PC I = PCNEW 
NNN=NNN+1 

IFINNN.GE.5) GO TO 18 
GO TO 9 
13 CONTINUE 
CALL OUTPUT 
WR IT E ( 6 * 30 ) PC IG 
OELT AT=OEL TT 
MOIS=DISM 
SU*AB=SUMABi 
PONOZ -PONO Z I 
RHEAD=R E EAD I 
RNUZ =RNOZ I 
PH LAD-P HEAD I 
MNUZ=MNOZ I 

IF( IPO.N6.2.AND. IP0.NE.3) GO TO 53 
NDUM = 1 
CALL OUTPUT 
NDUM® 0 
53 CONTINUE 
IPT = 0 
RETURN 
END 
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> 

> 

> 
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> 
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SUBROUTINE PL OT I T ( X , XHDR , KX, Y f YHOR » NY » T » THOR , NT ,NP, NPLOT .DUMMY, NO I 


******** ** ****** ******* ********************************************** 

* SUBROUTINE PLOTIT PLOTS TWO OE PENCE NT VARIABLES, Y AND T, * 

* VERSUS AN INDEPENDENT VARIABLE, X * 

* XHDR, YHPR, AND THDR ARE THE HEADINGS FOR THE X, Y, AND T * 

* AXES, RESPECTIVELY * 

* KX, NY, AND NT ARE THE NUMBER OF CHARACTERS IN THE X, Y, AND * 

* T AXES HEADINGS, RESPECTIVELY (MAX OF 32 IN EACH) * 

* NP IS THE NUMBER OF POINTS TO BE PLOTTED PLUS 2 * 


* 

* 

* 

* 

* 

v 

* 

* 


VALUES FOR NPLOT ARE 

1 FOR Y ONLY PLOTTED VERSUS X 

2 FOR Y AND T PLOTTED VERSUS X ON SAME AXES 

WITH INDIVIDUAL SCALES 

3 FOR Y AND T PLOTTED VERSUS X ON SAME AXES 

WITH SAME SCALES 

DUMMY IS THE HEADING FOR THE DOUBLE AXIS (NPLOT=3) 

NO IS THE NUMBER OF CHARACTERS IN DUMMY 


$*##***£*#***********#**« **#*#«#*«***£***** ^^tf#***********#*#******-* 


DIMENSION XHDR ( B ) , YHDR ( 8 ) , THDR ( 8 ) .DUMMY! 3) »X(NP)»Y(NP)»T(NP) 
CALL GSI2E (12.0,11.0,1121) 


NX =— K X 


NM=NP-1 
NN = N P- 2 

I F ( NPLQ T . EC • 1 ) GO TO 9 
CALL SCALE (T, A., JN ♦ 1 > 

9 CALL SC ALE ( X, B. ,NN , l ) 

CALL SCALE! Y, A. ,NN, 1 ) 

CALL PLQT( fc.25, 2. ,-3) 

IF (NPLOT. NE. 3 ) CALL AX IS l 0. , 0. , YHOR ,NY,4. ,180. , Y< NM ) , Y ( NP) ) 
IF (NPLOT. FQ. 3) CALL A X I S ( 0. , 0. ,0 UMM Y, NO ,4. , 1 80 . , Y l N K) , Y ( N P ) ) 
CALL AX IS(C.» 0., XHDR, NX, 8. ,90. »X(NM) , X(NP) ) 

IF(NPLOT.EQ.l) GO TO 12 


DO • U I = 1 , NN 

11 T( I )=-T( I I 

12 DO 13 I = 1, NN 

13 Y( I >=-Y (I) 

CALL LINE(Y»X»NN»1,L,1) 

CALL PLOT! G., 0., 3) 

IF(NPLOT.EQ.l) GO TO 24 

IF (NPLOT. EQ. 2 ) C A LL PLOT { 0 . . 5 , 2 ) 

IF (NPLOT .E0.2 ) CALL AX I S ( 0. . 5 • THDR , NT ,4. , 1 80 . ,T ( NM ) * T ( NP ) ) 
CALL LINE(T*X,NN,1,1,2) 

DO 2 5 I = 1 , NN 

25 T ( 1 )=-T( I ) 

24 DO 26 1=1, NN 

26 Y(I ) =— Y ( I ) 

IF(NPLOT.EQ.l ) GO TO 32 

CALL SYMBOL (-4.35, . 52,.1,1,0.»0) 
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A A A A A A 


AGE 32 


I 


TABLE IV-2. (Cont'd) 


CALL SYMBOH-4.2, .52, .1,2,0. ,0) 

CALL SYMBOL {- 4.3,. 65, .1, YHDR » 90. , NY) 
CALL SYMBOL (-4. 15, .65,. I, THOR, 90., NT) 
32 CALL PL0T(8.5,0.,-3) 

RETURN 

END 
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Figure IV-U. Flowchart For ignition subroutine. 
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Figure IV-5. Flowchart for plotting subroutine 






V. TEST CASES 


In this section two test cases of the extended computer program are 
presented. Complete input data for one test case and samples of the output 
data for both cases demonstrating the extended capabilities of the program 
are presented. 

Test Case 1 - Modified SUM 1 ✓ 


The SRM represented in this case is the same as SRM 1 of the original 
report (Reference 1, Figure V-l) except that the forward segment, a 
truncated star grain, was replaced by a wagon wheel grain defined in the 
table of input data (Table V-l, present report). Figure V-l demonstrates 
results of the ignition transients calculations for this SRM in the form of 
the actual plot obtained from the CalComp 663 digital plotter and Table V-2 
shows a sample of the associated numerical printout . Plotter output of 
additional results are given in Figures V-2 and V-3. 

Test Case 2 - SEM 2 

The SRM represented in this case is identical with SRM 2 of the original 
report (See Reference 1, Figure V-3 and Table V-4). However, to examine the 
new program capabilities , the performance of SRM 2 was calculated in two 
different ways : l) with the head end of the grain considered coni cal ( C0P=0) , 

and 2) with the head end considered hemispherical, the more accurate geometric 
representation, (C0P=3). The aft end in both cases was considered flat 
because the propellant downstream of the aft slot is not represented by 
equations but by tabular values . The results for the two cases are shown 
on -Figure V-4 and compared with static test results. It is notable that the 
two program results are significantly different. The hemispherical end 
representation gives a somewhat higher pressure level throughout equilibrium 
burning as would be expected from the nature of the geometric representa- 
tions . (in retrospect a higher value of TKETCH should have been chosen for 
the conical end approximation) 

Although the hemispherical end approximation yields results that are 
better than the conical end approximation for the first two-thirds of the 
trace, results for the last third are the opposite. Additional investiga- 
tion using the computer, program indicate that the departures from test results 
could possibly be the result of a generally high nozzle throat erosion rate 
(RADER) during the particular test than that used for the program computa- 
tion. 
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TABLE V-l. INPUT LATA FOB SBM 1 WITH WAGON-WHEEL IN PLACE OF 
TRUNCATED STAR GRAIN SEGMENT. 

OPTIONS 


[ 00= 1 
I WO = 1 
i P0= 1 

f'HjMPL T ( J J ) * 0» Of Ot Of Of Of Ot Of 
FRCPELLANI CHARACTERISTICS 
RHO* 0.001590 
A* 0.0712 
N= 0.250 
ALPHA= 0.0 
BETA* 0.0 
MU* 1.2700E-C7 
CSTAR* 5.166CE 03 

BASIC MOTOR DIMENSIONS 
L = 1243 . 

TAU = 50.47 
Oc- 1.A800E 02 
DII* 4 • 7COOE 01 

1 Ht: T A = 0.0 

A LF AN* 3.03406-01 
LTAP* 1 . 5000 1 02 
XT* 5.0000E 00 
70= O.C 

BASIC PERFORMANCE CONSTANTS 
OELTAY* 0.250 
XGUT* ICCO.CC 
C PCU T * 500.00 

2 E T AF * C.980 
T6 = 13C.C 
HB S 15CCOO. 

GAM* 1.143 
RADtR* C.0I35 

GRAIN CONFIGURATION 
INPUT- 3 
GRAIN* 3 
STAR* 3 
NT* C. 

ORDER* 1 

C0P= C C.P. GRAIN GEOMETRY 

EG* 154.940 

Cl* 54 . GOO 

CELDT* G.C 

$* 7. 

THETAG* C.C 
LGCI*' 1046. 50 
L G N I = C.C 
THETCN* C.75C35 
THETCH* 1.57075 

BASIC STAR GEOMETRY 

N 5= 1 . 

L 0 S I = 146.00 
N P = i 0 . 

PC* 77.47C 
FILL* 2.50 0 

NN= 0. 


0 , Of 0 f C f if o , Of 0 

WAGON WHEEL GEOMETRY 
TAUWW* 13.00 
L L* 30.00 
L 2= 1.00 

ALPHA 1 = 1.57060 
ALPHA 2= 1.57080 

HW* 10. GO 

IGNITION CCNS1ANTS 

KA* 5145. C 
KB* C.C1C 
UFS* 5C00.0 
CSIG* 49CC.C 
PM I G* I6CC.C 
T I 1= C.1C0 
T 1 2* 1.50 

RRIG* 172C.C 
DELTIG* C.CC5 
PBIG* 14.696 

INERT WEIGHT INPUTS 
PIPK* 1.5CC0E-03 
OTEMP* 3.CCCCE 01 
SIGMAP* 3.CCCCE-C2 
SIGMAS* 2.CCCCE-C2 
Nl* 3.CCCCE CC 
N2 * C.O 

SYCNOM* l.SOOOt C5 
CCC* 1.55506 02 
PSIC* 1.40Q0E 00 
DELC* 2.83CCE-C1 
LCC* i. 1025E 03 
NSEG* 5.C0CCE 00 
HCN* 6.2CCCE 01 
SYNNCM* 1.9C00E C5 
PSIS* 1.4CCCE 00 
PSIA* 2.0000E 00 
Kl* 2.C80CE-01 
K2= 5.250CE-02 

PS I INS* 2.COOOE 00 
OELINS* 4.6200E-C2 
KFH* 3.CC00E-03 
KEN* 1.2CCGE-02 
DLlNER* 3 . 3000 E- 02 
TAUL* 6.50Q0E-G2 
WA* C.G 

TABULAR VALUES 

VCI'T-0 

Y ABPK 
0.0 9,000 
9.67 9,000 

40.0 0 

160.0 0 

ALL OTHER A f s ARE 0. 
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TABLE V-2. SAMPLE PRINTOUT OF IGNITION COMPUTATIONS FOR MODIFIED SRM 1 


IGNITtR Sl/t CALCULATIONS 
A S IG = 52. <55 
WIGTOT= 404.62 
MIGAV* 10.785 


***********$** ***'***<** 1 M*< : <‘ 
**** IGNITICN TRANSIENT *<** 


T I Mh = 
R NOZ * 
P T AR = 
P A T M = 

PC IG= 


C.O 

c.o 

1 . 320 L £ 00 
1 .4696E G 1 
1.4696E Cl 


Y = O.C 


RHE AD= 

0.0 

PONOZ = 

1.46966 01 

PHEAC- 

1.4656E 

Cl 

MNOZ- 

0.0 

SUMA0= 

C.C 

SG= 

7.83316 

C2 

CFVAC = 

0.0 

F V AC = 

c.o 

F= 

0.0 



TIME- 0.025 
RNOZ = 1.41436-C1 

PTAR = 1.32016 00 

PA1N = 1.4696E 01 

PC IG= A .00006 C2 


Y= C.C 

RHL AD = 1.45556-01 

MNO Z « 5.21596-01 

C f V A C= 1.77756 00 


P ONOZ = 1.761CE 01 

SUMA B= 4.3135E 04 
F V A C- 5.21076 04 


PHE AC- 1.96626 01 
SG= 7.83316 C 2 
F= 0.0 



o GPHEfiD 

~ ^ PONOZ 



TIME (SECS) 


Figure V-l. Ignition transients for SRM 1 with wagon-wheel 
in place of truncated star grain segment. 



PRESSURE CPS I R) *1CT 

UO . 00 80.00 120.00 160.00 


G PHEflO 
> PQNQZ 



TIME (SECS) 


Figure V-2. Equilibrium burning and tail-off for SRM 1 with wagon-wheel 
in place of truncated star grain segment. 



PRESSURE CPS I R) *1 

.00 20.00 40.00 60.00 



Test data. Ref. 9 

Program, COP = 0 

— - — Program, COP = 3 


o 

V 


°cr.oo s', oo 


10.00 15.00 2o7ocT 25.00 

TIME (SECS) 


30.00 


35 . 00 


40.00 


Figure V-4. Comparison of test results for SRM 2 with computer results 
obtained using various methods of representing head end 
grain geometry . 


VI . CONCLUDING REMARKS 


The extended SRM design and analysis program presented in this report 
considerably improves the capability of the program by permitting calcula- 
tion of ignition transients and wagon-wheel type grain configurations . The 
plotting option offers an aid to the designer for rapid interpretation of 
the results. 

A number of changes in the program, notable among which are calcula- 
tion of initial and final gaseous chamber volume and more accurate geometric 
representation of the ends of circular-perforated grains should result in 
somewhat more accurate SRM performance predictions. Additional refinements 
are of course possible, but, as in the present case, the refinements will 
add to the complexity of the input preparation and the computer operating 
time. Before incorporating such changes, the degree of improvement in 
prediction capability anticipated should be evaluated in light of the basic 
objective of the present work of providing a simplified program and in light 
of the approximations inherent in the internal ballistics model used. 


- 64 - 



REFERENCES 


1 . Sforzini, R. H. , "Design and Performance Analysis of Solid-Propellant 
Rocket Motors Using a Simplified Computer Program," Final Report, 
prepared for NASA, Auburn University, October 1972. 

2. Sforzini, R. H. , and Fellows, H. L. , Jr., "Prediction of Ignition 
Transients in Solid— Propellant Rocket Motors," Journal of Spacecraft 
and Rockets , Vol. 7, No. 5, May 1970, pp. 626 - 628 . 

3. "Solid Rocket Motor Igniters," NASA Space Vehicle Design Criteria 
(Chemical Propulsion) NASA SP-8051, National Aeronautics and Space 
Administration, March 1971, p. 6 l. 


- 65 - 


NASA— MS FC 



